


AD A 0 1 5 1 4 0 


280089 


AFRPL-TR-75-36 


A COMPUTER PROGRAM FOR THE PREDICTION OF SOUD PROPELLANT 
ROCKET MOTOR PERFORMANCE, VOL. I 

FINAL REPORT 


Ultras /stems, Inc. 

2400 Michel son Drive 
Irvine, California 92664 


Authors: D. E. Coats 
J. N. Levine 
G. R. Nickerson 
T. J. Tyson 


N. S. Cohen 
D. P. Barry m 
C. F. Price 


JULY 19 7 5 



3 

f 


< 

T 


1 


APPROVED FOR PUBLIC RELEASE; 
DISTRIBUTION UNLIMITED 


D D 9- 

[ — \ r s r r- . 

Ui( stf s ’ s 


NATIONAL TECHNICAL 
INFORMATION SR Via 


1 « m| HI. VA. TDM 


AIR FORCE ROCKET PROPULSION LABOnAIORY 
DIRECTOR OF SCIENCE AND TECHNOLOGY 
AIR FORCE SYSTEMS COMMAND 
EDWARDS, CALIFORNIA 93523 


a. 




NOTICES 


“When U.S. Government d ravdegs, specifications, or other data ere used 
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tho Government thereby Insure no responsibility nor any obligation whatsoever, 
end the fact that the Government may have formulated, furnished, or In any way 
supplied the said drawings, specifications or other data, is not to be regarded by 
Implication or otherwise, or In any manner licensing the holder or any other person 
cr corporation, or conveying any rights or permission to manufacture use, or sell 
any patented invention that may in any way be related thereto. " • 


FOREWORD 


This report was submitted by Ultras y steins, Inc., Environmental and Applied 
Sciences Division, 2400 Mlchelson Drive, Irvine, California 92664, under Contract 
No. F046 1 1-73-C-0038, Job Order No. 305909LZ with the Air Force Rocket Propul- 
sion Laboratory, Edwards, CA 93523. 

This report consists of three volumes. Volume I describes a computer 
program for the prediction of Solid Propellant Rocket Motor Performance. The com- 
puter program described herein will be referred to as the SPP program, and describes 
the engineering analysis which was used in developing this computer program and 
the results obtained to date. 


Volume II of thl3 report is a programming document of the computer program 
which was developed under this contract. It Includes a subroutine-by-subroutlne 
description of all of the elements of the SPP program. 

Volume III of this report is a Program User's Manual which describes the 
input necessary to execute the SPP computer program and the information required 
to Interpret the output. A sample case Is also Included. 

This report has been reviewed by the Information Office/DOZ and is re- 
leasable to the National Technical Information Service (NTIS). At NTIS It will be 
available to the general public, Including foreign nations. 


This report is unclassified and suitable for public release. 
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A flexible, modular, fully automated, solid rocket motor performance prediction 
program has been developed. The program, which has been given the acronym 
SPP is based on six pre-existing computer codes. These codes have been Inte- 
grated and modified, as required. To supplement the theory, where necessary, 
and to increase the flexibility of the pregram, a number of existing and newly 
developed semi-empL-ical correlat:. rt \-cre Incorporated into the program. The 
program has a general three-dimensional grain design capability, cowled to a 
one-dlmenslonal ballistics analysis. Tie deviations from Ideal performance are 
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computed as a series of Independent efficiencies. Hie program currently treats 
the following losses: two-dtmensianal/two-phase (coupled), nozzle erosion, 
kinetics, boundary layer, combustion efficiency, submergence. The program 
predicts average delivered perfonnar.ee, as well as mass flow, pressure, thrust, 
impulse, end specific impulse as functions of time and trajectory. 

In order to assess the validity of the SPP program, calculated results were com- 
pared to firing data for four different types of motors. While conclusive state- 
ments regarding the accuracy and range of validity of the SPP program cannot be 
made until additional verification efforts are conducted, the results of these four 
test cases were encouraging. These calculations also served to demonstrate 
the desirability of eliminating some of the present limitations of the program. 
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PREFACE 





"his Is Volume I of a three-port report which describes a computer 
program for the prediction of Solid Propellant Rocket Motor Performance. The 
computer program described herein will be referred to as the SPP program. 

Volume I of this report describes the engineering analysis which was 
used In developing this computer program and the results obtained to date. 

Volume II of this report Is a programming document of the computer 
program which was developed under this contract. It includes a subroutlne-by- 
subroutlr.e description of all of the elements of the SPP program. 

Volume III of this report Is a Program User’s Manual which describes 
the input necessary to execute the SPP computer program and the Information 
required to Interpret the output. A sample case is also Included. 
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ABSTRACT 

A flexible, modular, fully automated, solid rocket motor performance pre- 
diction program has been developed. The program, which has been given the 
acronym SPP is based on six pre-existing computer codes. These codes have 
been Integrated and modified, as required. To supplement the theory, where 
necessary, and to Increase the flexibility of the program, a number of existing 
and newly developed seml-emplrlcal correlations were Incorporated Into the pro- 
gram. The program has a general three-dimensional grain design capability, 
coupled to a cne-dlmenslonal ballistics analysis. The deviations from ideal 
performance are computed as a series of Independent efficiencies. The program 
currently treats the following losses: two-dlraenslonal/two-phase (coupled), 
nozzle erosion, kinetics, boundary layer, combustion efficiency, submergence. 
The program predicts average delivered performance, as well as mass flow, pres- 
sure, thrust. Impulse, end specific Impulse as functions of time and trajectory. 

The theoretical models and empirical correlations upon which the program 
Is based are described, together with the assumptions that are employed. The 
known limitations of the analysts are presented to aid the user In determining 
the range of applicability of the program. 

In order to assess the validity of the SPP program, calculated results were 
compared to firing data for four different types of motors. While conclusive state- 
ments regarding the accuracy and range of validity of the SPP program cannot be 
made until additional verification efforts are conducted, the results of these four 
test cases were encouraging. Thtse calculations also served to demonstrate 
the desirability of eliminating some of the present limitations of the program. 
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1 . INTRODUCTION 


The physical phenomena which govern the overall efficiency of a solid rocket 
motor are quite numerous and complex. For this reason the prediction of solid 
rocket motor delivered performance has, historically, been empirically oriented. 
Total reliance on correlations of firing data Is undesirable, however, as their 
ranges of applicability cannot usually be clearly defined, and a priori performance 
predictions for new motor designs or propellants cannot be made with confidence. 

In order to achieve the desired generality, performance predictions should, where 
possible, be based on physically realistic models of the controlling phenomena. 

Over the past ten, to fifteen, years Individual facets of the performance 
prediction problem In solid rocket motors have been the subject of considerable 
study. Reference i contains a survey of the state-of-the-art up to 1971. Some 
of the more recent developments In this field (Including parts of the present Inves- 
tigation) are discussed In Reference 2. Terse description of rocket industry 
analytical performance prediction capabilities as of July 1974 may be found in 
Reference 3. Some recent analytical and experimental work is described In Ref- 
erences 68 and 69. 

When this effort was Initiated, analytical models for most of the physical 
processes governing motor performance had been formulated, and computer programs 
based on these models were In existence. Unfortunately, these computer programs 
were developed Independently; each program addressing a particular aspect of the 
overall oroblem. Few organizations possessed a large enough operational subset 
of these programs to make analytical performance predictions. This situation pro- 
vided the primary motivation for the present effort. 

It was believed that the existing technology was capable of yielding adequate 
performance predictions, but It was felt that a systematic and efficient procedure 
for Implementing this technology was required In order to foster Its use on a in- 
dustry-wide basis. Thus, the primary goal of this program was to develop an 
automated, computerized, methodology, based on existing methods, which would 
be capable of predicting delivered specific Impulse to within +2% and delivered 
thrust and total impulse to within +5%. The method of approach utilized to achieve 
this goal is outlined in Section Since we can expect continuing advances in 
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the state-of-the-art, the Solid Rocket Performance Program (SPP) was modularly 
structured so as to be able to easily accommodate them. 

The theoretical techniques employed In the SPP program are succintly de- 
scribed In Section 3. This report Is not meant to be a detailed dissertation on 
solid rocket performance prediction methods. It's purpose Is to Indicate the methods 
selected to meet the objectives of this study, the reasons for their selection, their 
primary assumptions; and limitations. Detailed descriptions of the theoretical 
models and solution techniques may be found In the pre-existing source material 
referenced throughout the text. However, In cases where existing methods were 
significantly modified, or extended, details of the changes are provided herein. 

Pre-existing and newly developed performance loss correlations are presented 
In Section 4. Certain of these correlations, e.g. combustion efficiency, are re- 
quired since theoretical models for all of the pertinent phenomena have yet to be 
developed; the others have been Incorporated Into the program for the sake of In- 
creased flexibility. 

No performance prediction methodology would be complete without considering 
such effects as bum rate, nozzle erosion and particle size. The manner In which 
these, and other Items, were addressed Is described In Section 5. 

The limited comparisons of analytical predictions with test data that have 
been carried out to date are described In Section 6. Conclusions drawn from this 
study are presented In Section 7. 1. Section 7.2 contains a list of the more signi- 
ficant limitations of the SPP program and our recommendations for future efforts 
directed tcvards alleviating them. This Is followed by Section 8, which contains 
the references cited In the text. 

As noted In the Preface this report Is Volume 1 of a, three volume set. Since 
this volume and Volume 3 (User's Manual) will often be used by different parties 
some of the descriptive material contained herein has also been Incorporated Into 
the User's Manual. 
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In accordence with die last two of the listed program characteristics the 
Solid Rocket Performance Program (SPP), was modularly structured; each module 
(with the exceotion of the control module) having been based on a pre-existing 
computer program. The following programs form the foundation of the SPP program: 

• One Dimensional Equilibrium Program (ODE)^ 

• One Dimensional Kinetics Program {ODKJ^ 

• Two Dimensional Two Phase Program (TD2P)® 

( 7 ) 

• Turbulent Boundary Layer Program (TEL) ' 

( 8 ) 

• Three Dimensional Grain Design Program 

• Motor Ballistics Program ® 

Each of the first four of these programs Is a separate module In Jie TPP pro- 
gram. while the last two programs have been combined In a single .nad'iic. These 
p rogra ms have been combined with a control module which permits th .an to be run 
together automatically. The six modules that comprise the SPP progra-i are sche- 
matically shown In Figure 2-1. The primary function, or functions, of each 
module Is also Indicated. All of the programs which were adopted for ^se In the 
SPP program have been modified to some degree to better suit the present purpose. 
These modifications are discussed In Section 3, as part of the module by module 
description of the program. 

The SPP program has also been designed to provide the user of the code 
flexibility In selecting which of the modules are to be exercised In any one com- 
puter run. The control module allows all jf the modules and programs to be ex- 
ercised sequentially In a single run. Often, however, this Is not the most efficient 
approach. Therefore, the program also allows the performance related modules to 
be exercised Individually, or In any "aer determined combination. A complete per- 
formance prediction, utilizing all of the program modules can be carried out in 
parts as efficiently (from a computer time standpoint) as If only a single run were 
made. This Is made possible by the nature of die linkage structure built Into the 
program. All of the data whit* needs to be transferred from one module to another 
is written out on an external logical unit (In addition to being automatically trans- 
ferred Internally). At the completion of a run, this linkage data can then be punch- 
ed, or stored In some other permanent manner. This data Is then available for use 
in subsequent runs, thereby eilmi-v >>e need for redundantly exercising a module 
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This type of linkage has two major advantages. The first is an obvious 
savings In computer run time for parametric studies in which the output from one or 
more modules is unchanged from case to case. The second is the ability to use 
other analyses Of appropriate) by 'fooling' the SFP code into accepting them as 
an integral part of its own computational procedure. This allows the user to sub- 
stitute die output from any other code or analysis into the SPP prefonnance predic- 
tion procedure, provided It is m Je to conform to die proper format. 

In conclusion, the SPP program may be described as a modular, flexible, 
code which combines theory and empiricism to provide accurate predictions of 
solid rocket motor delivered performance. The analytical techniques and correla- 
tions currently utilized in die SPP are summarized in Table 2-1. In many cases the 
program allows the user (an option) to directly input a performance related quantity. 
In these cases (indicated In die table) the pertinent analytical and/or empirical 
calculations are bypassed, in favor of the input quantity. 

2.2 Definition of tosses 

The basic method of approach used in the development of the SPP computer 
program is to calculate the maximum or ideal ‘theoretical’ performance of a rocket 
motor and to subtract from that maximum performance the losses or deviations from 
the "theoretical" performance which are known to occur in a real rocket motor. 

In Implementing such an approach It 13 often assumed that the various losses can 
be treated independently. As noted In Section 2.3, all of the losses should not be 
treated independently. There are some Interactions strong enough to warrant 
consideration. With these exceptions in mind, the concept of treating the losses 
Independently appears to be valid, considering the success of this approach, and 
others of a similar nature. 

A standardized performance nomenclature does not exist. Therefore, for the 
sake of clarity, the definitions used herein are listed below. 

• Performance - unless otherwise specified, performance 

will be in terms of specific Impulse I , 
delivered to a vacuum. sp 

•Theoretical I - the maximum possible which can be de- 

v ltvered at the initial area ratio by a propel- 

lant at a given chamber pressure and total 
enthalpy. Referred to as I, 

sp th 
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SUMMARY OF SPP PROGRAM CAPABILITIES 








• Loss 

• Two phase flow loss 

• Two dimensional or 
divergence loss 

•Finite Rate Kinetics 
loss 

•Boundary Layer Loss 

•Erosion Loss 

•Submergence Loss 

•Combustion Efficiency 

•Delivered Performance 


the decrement from the theoretical l £p which 
can be attributed to a physical phenomena 
not Included In the calculation of I.. . 

sp th 

the decrement In performance due to finite 
velocity and temperature differences be- 
tween the gas and condensed phase. 

the decrement In performance which Is due 
to the momentum of the rocket exhaust not 
being totally aligned with the axis of the 
motor. 

the decrement in performance which is due 
to Incomplete transfer of latent heat to sen- 
sible heat caused by the finite time required 
for gas phase chemical reactions to occur. 

the decrement In performance due to viscous 
forces adjacent to the nozzle wall and heat 
transfer to the nozzle wall. 

the decrement In performance associated 
with erosion Induced changes In the nozzle 
exit area ratio 

the decrement In performance that can only 
be attributed to the fact that a given nozzle 
is submerged. 

as used here, the term combustion efficiency 
will refer only to the degradation of I due 
to a departure of chamber total temperature 
from the theoretically calculated total temper- 
ature. 

the calculated performance of the rocket 
motor which Includes all of the losses con- 
sidered herein. 


2.3 Calculation of Delivered Performance 

Within fte framework already discussed (l.e., separation of Individual 

losses), there are three basic methods which can be used to calculate delivered 

I . The first method Is to calculate discrete decrements In performance, aI 
sp sp 

and to subtract these losses from the theoretical specific Impulse, I . That 
is Ptn 


sp D sp th 


- Tj Al, 


sp. 


loss 


( 2 - 1 / 
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The second method is to calculate a series of efficiencies, n's, to be applied 
multipllcatively, l.e. 


^Pp “ ^ " »>10S3 • (2_2) 

Both of the above methods are the same through first order, thus. If the losses 
are small they yield the same result. 

The third method, selected by the J ANN AT Committee on Performance Stan- 
dardization for liquid rocket engines. Is a combination of the first two methods. 
That Is, all of the losses except the boundary layer loss are treated as efficiencies 
while the boundary layer loss Is treated as a decrement. Hence, 

Isp D “ Isp th ” 1,1055 ' 4Iap BL <2 ' 3) 

This latter method of approach was selected for Incorporation In the SPP code for 
two reasons. The first reason being that there was no compelling reason to de- 
viate from what has been established as a “standardized" method, and secondly, 
that any other combination of the above methods yields substantially the same end 
result. 

Unless a grain design and ballistics solution Is obtained, the only other 
performance related quantity that Is calculated Is the specific Impulse delivered 

to average ambient pressure, I , defined as follows: 

sp D 

a 

i =: -L. 

sp D sp D P 
a c 

where: 

P a average ambient pressure (pst) 

Q 

P c = average chamber pressure (psl) 
r = average nozzle expansion ratio 
CJ, = theoretical C* (ft/sec) 

C D = nozzle discharge coefficient (from (TD2P) 


^ C ‘th 

-egr- 


(2-4) 


4 



'J 

3 


i 
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If a grain design and ballistics solution is obtained delivered thrust (vacuum 

and ambient), delivered specific Impulse (vacuum and ambient) and mass flow are 

calculated as functions of time. The I values are calculated by applying equa- 

sp 

ting (2-3) and (2-4) at successive time Increments, Delivered thrust values are 
then computed as follows*, 

F D (t) = m(t) I (t) (2-5) 

u sp D 

F D (t)= F D (t) - P a (t) A e (2-6) 

a 

The total mass flow, m^, total impulse 1^, time averaged specific Impulse, , 
and time averaged mass flow rate, m are also calculated in the indicated manneP, 


m T = 

J l bum m(t) dt 

0 

(2-7) 

Ij = 

fHjum F_ (t) dt 
o a 

(2-8) 

T *p d = 

V m T 

(2— S) 

5 

m '/ t bum 

(2-10) 


where tj Jum denotes the total motor bum time. 

2.4 loss Interactions and Reference Performance Values 

The concept of defining a theoretical I and various individual losses Is 

sp 

a convenient artifice that allows us to predict delivered preformance. In addition 
to being artificial, this concept of treating individual losses is somewhat sub- 
jective. The number of “losses" that are defined, and what they are called, 
can, and does, vary from one performance prediction methodology to another. 

There is a constraint on this subjectivism, however, which unfortunately is some- 
times violated. One way of expressing this constraint is as follows: any number 
of individual losses may be defined, as long as each physical source of non-idea! 
performance is accounted for once, and only once. 

•If a grain design and ballistics solution is obtained, but a two phase flow solution 
is not, then thrust is calculated somewhat differently. See Section 3.3.2 for the 
related details. 
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2.4.1 Loss Interactions 

In defining a set of losses one normally tries to arrange for each loss to 
be primarily determined by a single physical phenomenon or flow characteristic. 

Since the phenomena controlling motor performance are, in reality, a nonlinear 
mutually interacting set a one to one correspondence between phenomena and 
losses can be achieved only in a relative sense. Each of the "real", non-ideal, 
processes Influences each of the losses, but to varying degrees. In order to satisfy 
the aforementioned constraint all of the significant Influences of each phenomena 
should be accounted for, but only once. The losses listed In Section 2.2 were 
defined such that each Is primarily the result of a single physical phenomenon, 
or effect. For example, the kinetics loss Is chiefly a result of finite rate kinetic 
effects. Table 2-2 (from Ref. 10) represents an attempt to quantify the degree 
to which a given phenomenon Influences the losses that are not directly related 
to It. Most, but not all of the losses and phenomena are listed in this table. 

The table should be read down and to the left. The number ihat appears In each 
square (except the diagonals, which represent the Interaction of a loss with Its 
chief determinant) Indicates the relative degree of Interaction between the Indi- 
cated loss and phenomenon. The key to Interpreting this number is given In the 
legend below the table. The Importance of the various Interactions Is thus char- 
acterized In terms of their potential for Influencing predicted delivered performance. 

In developing the present methodology It was decided to directly Incorporate 
only those Interactions which were felt to be of primary Importance. It can be seen 
from T.-ble 2-2 that there are two such Interactions. These are: the effect of non 
one-dlmenslonal flow on the boundary layer loss; the effect of multiphase 
flow on the divergence loss. 3 

Two dimensional effects have an Important Impact on the nozzle boundary i 

layer because of the sensitivity of the boundary layer to free stream conditions. j 

The wall streamline In a typical nozzle differs substantially from the one dimen- 
sional analysts streamline for the nozzle, particularly In the nozzle throat region 
where the heat transfer rates are the highest. 

The effect of multiphase flow on divergence loss Is Important because of die 
manner In which heat Is transfered from the particles to the gas. Because of their 
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inertia, particle concentrations will be greater near the nozzle axis than the nozzle 
wall. Consequently, there will be greater heat addition to the axis streamline 
than to the wall streamline. This Is demonstrated by the fact that two dimensional 
two phase flow calculations for a perfect gas show a greater maximum velocity 
for the axis streamline than for the wall streamline. Thus, two phase flow effects 
cause a reduction In the percentage of misaligned thrust (lower divergence loss). 

It should be pointed out that the estimates given In Table 2-2 are based on 
typical motor configurations and operating conditions. The Importance of the various 
Interactions should be reexamined In extreme cases. For Instance, In a motor 
operating at very low pressure, with a metallzed propellant, the coupling between 
Incomplete combustion and the kinetics loss could become significant. 

2.4,2 Reference Performance Values 

Some of the Individual loss calculations directly yield a performance decre- 
ment or an efficiency, e.g. , the boundary layer loss. The other performance loss 

calculations yield I values, e.g. kinetics loss, 2-D, 2-phase loss. It Is not 
sp 

desirable to directly utilize any of these I values In the calculation of delivered 

performance. Rather, these calculated I values should be related to a reference 

sp 

value of I ; thereby producing either an efficiency, or a performance decrement.* 

Properly defined reference values are one of the keys to successful implemen- 
tation of a performance prediction methodology based on separation of losses. 

One must make sure that the reference computation Includes all of the same physical 
phenoir :oa and assumptions as the related loss calculation, except for the phen- 
omenon, or phenomena, responsible for the loss, as defined. Ideally, this reference 
calculation should be made by the same computer code which Is selected to make 
the loss calculation, since this would minimize deviations due to numerical tech- 
niques, truncation and roundoff errors. However, such an approach Is not always 
practicable. For example, both the equations for the two phase flow and chemical 
kinetics losses become "stiff” In the equilibrium limit and result In unacceptably 
long and difficult calculatlonal procedures. Even 11 sophisticated numerical tech- 


*ln the present methodology the use of efficiencies is preferred. It Is felt that 
efficiencies better compensate for the omission of things like real gas effects. 
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% ' nlques are used to calculate the "near" equilibrium case, extrapolation to the full 

| j equilibrium solution introduces errors which can negate the advantage gained by 

using the same computer code. 

In developing the reference performance calculations used in the present 
methodology the following criteria were used: 

The conceptual differences between the reference and loss calculations 
should Include only those mechanisms which cause the loss. 

Computational efficiency should be considered, without sacrificing accuracy 
and consistency. 

The modifications to existing computer codes should be Kept to a minimum 
to minimize program development costs. 

The adopted procedures are described in Section 3 . 



2-12 


3 


ANALYTICAL METHODS 
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The *heoreticai methods employed In the SPP program were all adapted from 
existing computer programs. Thus, the basic theories ard computational tech- 
niques have all been previously documented. As a result of the number of analyses 
contained In the present methodology no attempt was made to combine these pre- 
existing documents Into a single entity. The purpose of this section Is to indicate 
the methods selected, the reasons for their selection, their primary assumptions 
and limitations. Wherever applicable, details of current modifications and/or 
extensions to the pre-existing analyses are Included. This section Is structured 
along the lines of the program, i.e. the methods are discussed In relation to the 
module In which they are utilized. The Interfaces between the modules are also 
described. While It does not include analytical methods, in the usual sense of 
the word, the Master Control Module Is described first, since It Is a key element 
in the overall program structure . 


3.1 


Master Control Module 


The master control module (MC) controls the execution of the SPP computer 
code. The MC module selects, via user Input, which of the five basic calculation 
modules aie to be executed, and whether calculated, input, or empirical, losses 
are to be selected for the calculation of delivered specific Impulse, thrust, and 
total Impulse. This module also controls all of the internal communication bet- 
ween modules. 


In addition to its control functions the MC module stores selected results 
from each of t*-e performance loss modules and then uses them to calculate delivered 
performance. Table 3-1 shows the variables used to calculate delivered perfor- 
mance, and Indicates the module In which each Is calculated. 

1 The delivered performance calculation does not use all of the variables 11st- 

- ed In Table 3-1 In their given form. 
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Where necessary the MC mocule combines, ratios, averages, etc., these 
variables, in order to generate the desired quantities. The empirical performance 
loss estimates are also computed in the MC module. 
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Table 3-1. Variables Used In Delivered Performance Calculation 


\ 1 " ■■ - 

Module 

Variable 

Theoretical Performance 

Sp th sp REQ 1 

Grain Design and Ballistics 

P c , m, r*, L*, n c * as functions of t. 

Kinetics 

Is pkin at '»• 

TWo Dimensional TWo Phase Flow 

I , as functions of c and C_. i 

sp TD2P ***** u ‘ 

Turbulent Boundary Layer 

£1 

s Ptbl 



The equations used to calculate delivered performance have been given In 
Section 2.3. The equation fcr vacuum delivered specific Impulse, previously written 
in symbolic form (equation 2-3), is written here In the manner used In the program. 

'sp D = ! sp th ’taN ' t, TD 2P * ’’C. * 'ISUB * ' ! sp m (3_1) 


The theoretical specific Impulse, I 


Is losed on the Initial (t=0) nozzle 


exit area ratio. 


The loss due to finite rate chemical kinetics Is treated as an efficiency 
factor, The kinetic efficiency is calculated as the ratio of an I calcula- 

ted by the CDK module, to an I calculated by a special option In the ODE module. 

sp 

This special ODE option, referred to as the "Restricted Equilibrium" option causes 
the ODE program to compute an I based on the more restrictive physical assump- 

Sp 

tlons employed In the ODK calculation. This procedure, which Is further described 
In Section 3.4, Is r.ecessary if certain losses are not to be counted twice. 

If a motor experiences throat erosion during the course of a firing the nozzle 
exit area ratio, and hence, the specific impulse, varies with time. The performance 
decrement associated with this phenomena Is usually called the erosion loss. 

While this loss can be estimated based on one-dlmenslonal calculations. It was 
felt that a two-dimensional estimate of this effect would be superior. Thus, for 
eroding nozzles, ls calculated as follows: 
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* P TD2F 
Cp ID2P *«!* 


The quantity t^q 2 p ls time averaged two dimensional two phase efficiency. 
The TD2P module calculates an efficiency as a function of the area ratio t . Area 
ratio js a function of time Is established from the Input, or internally calculated, 
nozzle erosion rate. This allows nxD2 P to *■** calculated as 



As defined In equation (3-1) Vjjy 2 p Includes the 'erosion* loss. In addition 

to the coupled two-phase and divergence losses. The erosion loss, as defined 

herein. Is equal to the ratio of the 2-D, 2-phase I at the time averaged expan- 
— ®P 

slon ratio, 1 ( t ), to the specific Impulse evaluated at the Initial expansion 

“ P TD?P 

ratio, I ( (,). (Future modifications to the program should break the erosion 

sp TD2P 1 

less out as a separate efficiency.) 


The combustion efficiency, jjqj. Is obtained from the average empirical c* 
efficiency as follows: 


Vet = C 'D V c * 


(3-4) 


tjce = 1.0, If the previous relation yields a value 
greater than unity 

The c* t. clency Is a correlation cf experimentally measured c* values (see Sec- 
tion 4.1) and, therefore. Includes all of die physical effects which determine the 
nozzle mass flow rate. One of these effects, the combined two dimensional, two 
phase, effect on moss flow, is also Included In the calculation Gf Vye2P‘ T *’ e 
definition of jj C£ In equation (3-4) prevents the mass flow effect related to the nozzle 
discharge coefficient, C^, from being counted twice. 

It is well known that the discharge coefficient for single phase, two dimen- 
sional, flow Is less than unity. However, It is also well established that In non- 
equilibrium two phase flew more mass flows through a nozzle at a given chamber 
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pressure than would If the particles were In velocity and thermal equilibrium. By 
Itself this effect would produce C Q values greater Qian unity. For most metallzed 
(olid propellant motors the net result of these two opposing effects is a discharge 
coefficient greater than unity. Since ij c * is empirically obtained, and subject to 
scatter, the second part cf equation (3-4) Is added to prevent meaningless values 
of tiggfroa being calculated (i.e. t?cE>U- 

If the combustion efficiency were 100% {complete combustion) the c* efficiency 
of the motor, rj c *. would he given by 

n c . = &-S) 

Thus, in such cases, tj c * would usually be less than unity for metallzed propellants, 
and would always be greater than unity for non-metalized propellants.* 

The submergence efficiency, r) sus , is based on an empirical correlation and 
is a function of the throat radius, mass fraction of condensed phase, length of sub- 
mergence to length of the internal motor, and chamber pressure. This correlation 
is discussed in more detail in Section 4.2.6. 

The decrement in performance due to the formation of a turbulent boundary 

layer in the motor nozzle is transmitted directly to the MC module and is reported 

as a function of area ratio. The decrement in I Is taken to be the value corres- 

sp 

ponding to the initial nozzle exit area ratio. This decrement is assumed to be in- 
dependent of erosion induced area ratio changes, and throat roughness. 


Delivered performance (I only) is calculated somewhat differently when 
5 P 

based completely on the performance loss correlations presented in Section 4. 
Instead of equation {3»1), the following relation Is employed, 

= ’’KIN ^ 2 D ^2t Vcz ’’SUB " 4l sp & (3 ~ 6 


•Most solid rocket motors have normalized radii of curvature w 2 and, in the ab- 
sence of particles Cq tends to be .99 f Cp < 1 , or, in other words, with perfect 
combustion ri c * would be bounced as follows, 1 <i?c*S 1.01. This small effect 
tends to be masked by measurement and data reduction errors, and many times 
test engineers are reluctant to report values cf Jj c * greater than unity due to tho 
mistaken assumption that such values are physically Impossible to achieve. 
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In the simplified performance prediction methodology the two dimensional, 
tr 2D , and two phase, tj 2? , efficiencies are treated separately. The combustion 
efficiency, T!q£, Is simply equated to the empirical c* efficiency, rt c *! and the 
so-called 'erosion loss* Is not directly accounted for. If an estimate of the ero- 
sion Ion Is desired, and a TD2P solution haa not been obtained, one can compute 
the ratio («)/!.„ (rit^O)}, where I gp {< (t=0)} is the specific impulse at 


*P, 


th 


•P, 


th 


th 


the Initial nozzle expansion ratio, usually denoted simply by I 


sp, 


9i 


In order to 


sp. 


th 


compare the delivered specific Impulses obtained theoretically (equation (3-1)), 
and empirically (equation (3-6)) an estimate of the erosion loss must be addend erf 

to (3-6), or equivalently, I„ must be replaced by 1 (7). 

£p th sp th 

the Master Control Module also calculates the other performance related 

quantities of Interest I , F n ft), I*, etc. , using equations (2-4) to (2-10). 
sp D °a 

a 

3.2 Theoretical Performance Module 
3.2.1 Approach 

For a specified propellant combination the maximum specific impulse. I, 
which can possible be delivered can be predicted using any of a number of 
standard computer codes developed for this purpose. A survey conducted by the 
JANNAF Performance Standardization Working Group found sixteen such computer 
programs, each of which represented an independent development. All of these 
programs are. In principle, capable of calculating the quantities of Interest, i.e. 
shifting and frozen I gp and C p , exhaust species etc. However, many of these 
programs. In practice, cannot adequately handle the solid phases required In 
solid rocket motor applications. Of the remaining applicable programs, two stand 
out by virtue of their excellent characteristics. These programs are: 

The ODE Computer Program developed by NASAARC and 
described in NASA SP 273, Reference 4. 

The Theoretical Igp Program developed by the AFRPL and 
described In Reference? 10 and 11. 

with either computer program area ratio expansion calculations 


li 
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yielding theoretical performance can be carried out for any rocket engine of 
current or projected use provided basic thermodynamic data exists for the pro- 
pellant. Each computer program currently contains a vast library of thermodynamic 
data. A comparison of these two computer programs as applied to the proposed 
effort is presented in Hie paragraphs below. 

the OOE Program 

This computer pregram has been under continuous development at NASA/IRC 
since the mid 19S0's. In 1968 it was adapted by the ICRPG (currently JANNAF) 
Performance Standardization Working Group as the standard computer program for 
theoretical performance computation. A newer extensively revised version of ODE 
is now available. This version Is documented In complete detail in Reference 4. 

Older versions of ODE were unable to successfully compute equilibrium 
compositions fer aluminized propellants. This difficulty lead directly to the de- 
velopment of certain of the sixteen chemical equilibrium codes previously men- 
tioned. The more recent versions . the ODE program, however, are all capable 
of handling multiple solid phase products. For example, the SP273 version of 
ODE can successfully compute compositions with large amounts of solid phase 
carbon; a system with a history of difficulties. 

The ODE program uses a free energy minimisation technique, along with the 
following assumptions, to calculate the vacuum I and related theoretical per- 
formance parameters. 

One dimensional flow. 

Chemical equilibrium between Species and phases. 

Thermal and dynamic (velocity) equilibrium 

between condensed and gaseous phases. 

Adiabatic, Isertroplc flow. 

The ODE program requires curve fit data for enthalpy vs. temperature. Seven 
coefficients faj, a 2 , 2j, a^, a 5 , ag, a^) are required such that: 


CP 

p T 


a l + ®2 T + ®3 P + a i P + a S T * 


3-6 


aj + a 2 T/2 + a 3 T»/3 + a^/A + a g TVS + a g /T 



9j InT + a 2 T + agT 3 /2 + a 4 T 3 /3 + a s T*/4 + a y 


Two set* of coefficients are required per species such that: 

300 °K sis 1000°K for set 1 

1000 °K s T s 6000°K for set 2 

The fits must be continuous across the Juncture at 1000°K. NASA SP273 
contains coefficients for 421 chemical species. A NASA/LRC program exists for 
the purpose of calculating thermodynamic data and outputting the above curve 
fit data. This program, called PAC, Is described In Reference 23. A new all 
FORTRAN IV version of PAC Is now available and has been found to be highly sat- 
isfactory. 

The curve fit data Is highly compact and computationally efficient but suffers 
from two major draw backs, which are: 

• certain data cannot be accurately fit with only seven coefficients, 
especially over the range 300°K to 1000°K. 

• curve fit data Is awkward to produce since It requires use of the 
PAC program followed by careful comparison of the resulting fit 
to the original data. 

NASA SP273 also contains elemental composition and enthalpy data (l.e., 
reactant cards) for 62 commonly used propellants. 

Options available with the ODE program are: 

Point Options: 

• pressure and temperature 

• enthalpy and pressure 

• entropy and pressure 

• density or volume and enthalpy 


Problem Options: 

• rocket option (l.e. chamber, throat, assigned area ratios, and/or 
pressures) 
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AfRTI. Theoretical I Program 

This program also uses a free energy minimization technique to calculate 

theoretical performance, and employs the same assumptions listed for the ODE 

program. The Aran, program contains several options not found In ODE, e.g. 

Air Augmented I , Delta Velocity, Weight Optimization; however, these options 
sp 

are not of Interest for the present purposes. 

The Aran, program requires curve fit data for specific heat vs. temperature. 
Ten coefficients (Zj, z^, Zj, z^, z g , z g , z ? , z g , Zg, Zj Q ) are required such that 

C P = Z 1 + z 2 T ‘ 1 + z 3 T' 2 + z 4 T' 3 + z 5 T‘ 4 

for 298°KsTsl200°K 
and 

Op » z 6 + z 7 T + a 8 T*. z 9 T s + z 10 T* 

for 1200°K sir 6000°K. 

A curve fit program and an extensive library of fit thermodynamic data Is 
maintained at AfRPL. 

Either one of these two programs would be very satisfactory for the purpose 
of computing theoretical performance . Computer run times for both programs are 
short and are In no way a limiting factor. Computer core storage requirements are 
also modest (approximately 32K words) fci both programs. Both programs feature 
excellent input methods. The ODE program was chosen for incorporation into 
the performance prediction methodology. The deciding factors were as follows: 

1. The AFRPi. program lacked the extensive documentation provided with 
the ODE program. 

2. The ODE program has been adopted as a standard by JANNAF and Is widely 
distributed. 

3. The ODE program was already Integrated Into the Or.e Dimensional Kinetics 
Program (ODlQ , which was adopted for the kinetics loss calculation. Re- 
placing ODE with another program would have required a significant effort 
to link the new program to ODK 
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3,2.2 Modifications TO The ODE Program 


For the present purposes the ODE program was modified to calculate the 
transport properties required by the Grain Design and Ballistics Module, 2-D, 2 
Phase Module and Boundary Layer Module. In addition, the program has been 
modified to calculate the reference condition for the loss In performance due to 
finite rate chemical kinetics. This calculation Is referred to as the "restricted 
equilibrium’ solution and Is discussed In more detail In section 3.4. ODE has 
also been modified to perform an equilibrium calculation at the exact area ratio 
corresponding to solidification. To reduce the length of the program the subroutines 
not required for the performance calculation were deleted. Thus, only the so- 
called rocket problem option of the ODE program can be exercised In the present 
context. 

The transport properties calculated by the modified ODE program are the 
gas phase viscosity u , thermal conductivity, y , and Prandtl number, Pr. The 
viscosity and thermal conductivity of the individual gaseous species are calculated 
from formulas given In Reference 20: 

4.15822xl0' 8 , £ 7t 


(.45 + 1.32 




The required Leonard Jones parameters, Or^.f^) are Internally stored In the computer 
program for 206 species. (See Subroutine MUK, Volume II , for a list of the species.) 
The vlccoslty of the mixture Is calculated from Wilke's seml-emplrical formula^ 21 ) 


.-«f 
1=1 1 


“l U % *11 

1/1 
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where N Is the number of species, X^, the mole fraction of species 1, and Is 
defined by: 

, f M. l' 1/2 f U. 1/2 M, I/4I2 


, M w 1/2 M. 1/4 2 

*tf = p72 1 + ^ _ 1 + ( (i7 ) ( M j -5 (3 " 10) 

The thermal conductivity, *, Is based on the equation given by Mason and Saxena 
which Is a slight modification of Eucken*s relation. 


N N x. _1 

= S », 1 + 1.0652 -1- 

1=1 1 ]=1 x l 

1^1 


The Prandtl number Is simply given by 


A subroutine, 0MEGA, which performs a least squares fit for viscosity In the form 
(T/T*) w has also been added to ODE. This routine was necessary because 
the TDZP and TBL programs require viscosity to be In this form. The curve fits 
treat the chamber, throat and exit points; the throat values are treated exactly. 

3.3 Grain Design and Ballistics Module 

3.3.1 Grain Design Calculations 

3.3. 1.1 Approach 

It vms required that the SPP program be capable of calculating solid rocket 
motor delivered performance (l.e. , pressure and thrust versus time as well as 
specific Impulse and total Impulse) for a wide range of motor designs and grain 
designs. At a minimum, end-burning, circular port, star port and spherical geo- 
metries, and at least two other qeometrles representative of modem production 
motors (e.g., slotted tube and flnocyl), were to be considered. 

The problem with constructing separate specialized subroutines for each 
geometry Is that It Is both cumbersome and limited. A multiplicity «f Individual 
geometry subroutines would add greatly to program length. The resultant program 
would also be limited, since the list of potential grain geometries is exhaustive. 

In addition, the geometry calculations do not stand alone; they have to be Integra- 
ted with the ballistics calculations. v separate, specialized geometry cal- 




culatlon would have to be Individually Integrated with the ballistics calculations. 
Clearly, a general method, that accomodates a wide range of geometries. Is a 
more desirable means of fulfilling the objectives of this program. 

Several grain geometry computer programs of more or less general capability 
are available. Some are limited to the geometry problem, some Incorporate com- 
bustion and gasdynamtcs to provide an Integrated Interior ballistics capability. 
These candidate programs were reviewed In order to select the most appropriate. 
Selection criteria were principally geometric capability and adaptability to general 
Interior ballistics analysis. 

3. 3. 1.2 Review of Existing Computer Programs 

3. 3. 1.2.1 Aerojet Basic Grain Design and Interior Ballistics (564) 

This program solves the grain geometry and Internal ballistics problems 
In a sequential manner. The geometric solution Is two-dimensional, and the In- 
| temal ballistics uses a simple burning rate law and cne-dlmenslonal gasdynamtcs. 

Two-dimensional grain design can be very useful for symmetrical motors 
with sufficiently large length/diaraeter ratio that end effects become negligible. 

- However, at lower length/diameter, or with end surfaces burning, or with partial 

bum surface restrictions, even a nominal two-dimensional geometry would re- 
quire a three-dimensional treatment. Consequently, this program was not selected. 

3.3.1. '..2 Lockheed General Grain Design Program (470) 

This grain geometry program Is basically two-dimensional, but can accom- 
modate cylindrical grain designs In three dimensions. An Important feature of this 
program Is that It includes logic to account for nonuniform burning for ballistics 
adaptability. However, It was not selected because it is not fully capable of three- 
dimensional treatment. 

3. 3. 1.2. 3 Thiokol Generalized Grain Design Program (AGDA) 

The AGDA program is a three-dimensional grain geometry program, but Is 
not capable of handling Internal ballistics. The program uses the Elkton-Moore 
method which Is a vectorial approach to the geometric problem. It Is capable of 
handling a broad range of grain oo^igns, the complexity manifesting Itself In the 
effort required for Input. 
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This program was not selected because It does not Incorporate combustion 
and gasdynamlcs analysis, and does not contain logic for coupling such analysis 
to the geometry calculations. 

3. 3. 1.2. 4 Boeing Internal Ballistics Computer Program (D2-125286) 

The Boeing program Is a three-dimensional grain design program Incorpor- 
ating one-dlmenslonal gasdynamlcs and erosive burning. However, it was created 
for a special purpose and Is therefore deficient In some respects and excessive In 
others. It Is deficient In that It Is restricted to certain types of grain geometries. 
It Is excessive In die detail given to transients and high acceleration effects. 
Consequently, the program as presently constructed would not be best suited for 
the present purpose. 


3. 3.1. 2. S Hercules Grain Design and Internal Ballistics Evaluation Program (64101) 

This program consists of two parts, the Basic Grain Design Program and 
the Internal Ballistics Evaluation Subprogram. 

The Basic Grain Design Program examines a grain In three dimensions by 
analytical simulation of drafting techniques used In developing a grain design. 

It solves for geometric parameters used In Internal ballistics calculations, and 
also for mass properties used In flight simulation calculations. It was developed 
In order to handle the most complex geometries used In modem production motors, 
and therefore has the desired general capability. 


The Internal Ballistics Evaluation Subprogram uses simple one-dlmenslonal 
gas-dynamics and erosive burning. It uses the geometry data computed for each 
burning station In the Grain Design Program, and feeds back the ballistics data 
at that location. It therefore accommodates nonuniform burning, and contains 
logic for coupling grain design with Interior ballistics. Ignition and taUoff tran- 
sients are computed, but more for appearance than accuracy. 

This program was selected because of Its geometric generality and bal- 
listics adaptability. TWo conectlble deficiencies were detected: the method 
of computing the radial slots In segmented grains, and the enormous size of the 
program (largely due to subroutines not required for the present purpose). 
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3. 3. 1.3 The Selected Program 



The Hercules 64101 program Is described In detail in Reference 8. 'With 
all of its options, it consists of six links and sixty subroutines. Only a fraction 
of these subroutines, fifteen in all, were found to be necessary for the requirements 
of the present program. Most of these are contained within LINK 2 of the original 
program, which deals with the geometry calculations. LINK 4 (a pseudo ballistics 
option) and LINK S (plotting routines) were eliminated since they could not be ex- 
ercised when LINK 3 (the ballistic- - ubroutlnes) is used. The mass properties 
subroutines within LINK 2 were also eliminated. LINK 3 was modified to incor- 
porate the LPC ballistics analysis described in Section 3.3.2, and LINKS 0 and 
1 (main driver and common regions between LINKS 2 and 3) were modified in ac- 
cordance with the other insertions and deletions. These modifications resulted 
in an 80 percent reduction in the si 2 e of the program. Information regarding the 
operation of this modified program is presented in the Users' Manual portion of 
this report. Volume 3. 

Basically, the geometry calculations compute volumes, and changes in 
volumes. Other geometric quantities are derivable from these volumes. Thus, 
for example, the burr, area is equal to the rate of change of volume divided by 
the rate of change of web: 


m 

= ° P V 


(3-13) 

dm 

■ar 

= „ A. *£. 
°p"b dt 


(3-14) 


_ dV/dt 
dw/dt ~ 

4MJ.I 

fiw/6 1 'mean value 

(3-15) 


m = propellant mass 

m = propellant mass flow rate 

9 „ = propellant density 
P 

Ajj = propellant bum area 
r = propellant bum rate 
w = propellant web 
t = time 

V = propellant volume 


3-13 



I 


The changes In web over a time increment are calculated from the ballistics 
subroutine, and the corresponding changes in volume are calculated from the 
geometry subroutines. Tha computations Iterate on flow rate until a solution for 
pressure converges Bum area is output as a matter of information, but is not 
used as such In the computations. Volumes form the basis of the computations In 
order to best satisfy the constraint that the weight expended must equal the weight 
loaded. 

Inputs to the geometry calculations initially consist of the outer dimen- 
sions of the propellant grain, which is taken to be completely filled with propellant. 
Next, a series of inputs subtracts all initial void volumes representing the motor 
cavity. Such volumes may be represented by combinations of cylinders, cones, 
spheres, and prisms located by spatial coordinates; these volumes may Intersect 
each other, but must not Intrude into propellant. Supplementary parameters are 
available to define an Inhibited surface, a comer round, a segmented grain slot 
face, and the symmetry of the grain. The grain Is then divided Into an X-Y com- 
putational mesh. Web fractions are also defined. When the web burns past a 
specified fraction, information may be output or Internally transmitted between sub- 
routines. The finer the mesh, the greater the accuracy, the smoother the output 
and the longer the computer time. The method by which the changing volumes are 
calculated for each mesh segment Is presented in Ref . 8 and will not be repeated 
here. 

3. 3. 1.4 Program Modifications 

As far as the geometry calculations per se are concerned the only modifi- 
cation required was In the method of computing the volumes In the vicinity of 
opposed-face radial slots separating segmented grains. In the course of using 
the Lockheed 156-5 motor (a multiple segmented motor) as an early test case, 
there were pressure surges in the ballistics output. These were traced to exces- 
sive volumes being expended when the slots crossed the reference bum Increments. 
In this type of situation, the calculations were expending triangular wedges of 
propellant (formed by the end face and the port to a certain length) rather tnan 
orthogonal webs. The error was magnified by the multiplicity of the segments. 

This deficiency was corrected by providing Input logic to key the existence of a 
slot, and a single equation for the annular geometry of the slot to be used at 
that location. 
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A related problem was encountered with end-burning faces In general. 
Depending upon the selection of the mesh size and bum Intervals, It Is possible 
to have the same volume expended counted several times. This can be avoided 
by maintaining a certain relationship between these Inputs, as discussed In the 
User's Manual. No program modification was required. 

The aforementioned 80% size reduction from the original Hercules program 
was successfully accomplished. Care had to be exercised to assure that deletions 
would not destroy functionality. The retained subroutines are listed In the Pro- 
gram Description, Volume 2 of this report. 

The revised ballistics analysis Is discussed subsequently. 


3.3.2 Ballistics Calculations 


3.3.2. 1 Approach 

The Interior ballistics analysis solves the conservation equations applied 
to the rocket motor cavity to predict pressure and mass flow as a function of time. 

In its most rigorous form, a ballistics calculation solves the time dependent 
equations governing conservation of mass, momentum and energy. The simplest 
type of ballistics calculations, widely used In rocket motor design, employ a 
steady-state mass balance. The sophistication required Is determined bv the 
purpose of the analysis and the accuracy required. 

A comprehensive model of the energy exchange processes within a rocket 

motor was prepared by LPC for the prediction of duty cycles of controllable solid 
( 13 ) 

rocket motors . The sophistication was warranted by the emphasis on transient 
phenomena: predictions of Ignition delay times, ignition pressure spikes, oxtlngulsh- 
ment-relgnltion phenomena and transient performance Integrals were required. The 
program would be recommended for such purposes. For present purposes however, 
which emphasize steady-state performance prediction. Its size, complexity and run 
time would not be justified. 

Computer programs such as Hercules 64101, Aerojet 564 and Lockheed 
241 use a simplified one-dlmenslonal treatment. Although the Hercules 
64101 program was selected for Us geometric capabilities. It would not be re- 
commended for the computation cf Interior ballistics. It employs a three-stage 
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treatment with separate analyses for Ignition, steady-state and tailoff. The method 
of interfacing these stages is ill-defined. The transient stages do not employ an 
interna' flow analysis, and the steady-state stage uses a simplified one-dimensional 
flow solution with time dependent terms omitted. A better approach, to avoid ar- 
bitrary division of the motor operation into stages, would be to employ one set 
of equations. The Seeing program takes this approach, but like Reference 13 was 
developed fer a special purpose. 

The most recent interior ballistics analysis, completed by Lockheed for 

(9) 

the Air Force, resulted in the Lockheed 637 program* . This program couples 
a comprehensive ballistics analysis to a limited grain geometry capability. How- 
ever, the methodology is compatible for use with the Hercules geometry analysis 
for more general capability. Therefore, this analysis was selected as best-suited 
to the requirements of the present program. The effort required was to reconstruct 
LINK 3 of the original Hercules program, using the ballistics analysis from Ref. 9 
as the basts. The original designation for the ballistics subroutine In the Hercules 
program, MAIN3, was retained, 

3. 3. 2. 2 Flow Equations 

The equations for one-dlmenslonal flow used tr. the computer program are 
based on the following assumptions: 

Particles occupy negligible volume compared to the gas. 

Particle velocity and temperature are equal to the respective gas values. 

TV oartlcle-gas mixture behaves as an ■equivalent" perfect gas. 

Mass is generated from the propellant without an axial component of vel- 
ocity; although Invalid for end-burners , the error will be negligible be- 
cause of the large flow area downstream of the end face. 

Transient effects ure simply approximated. 

With these assumptions, the consecration equations are: 

P = e RT (3-16) 

iM + ifjA-*;, (3-17) 
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(3-18) 

h 2g 

• c p T r 

(3-19) 


P = pressure, llv'ft* 

p *» gas density, lb/if 

R « gas constant for the products, P/F 

T »•> gas temperature, °F 

t - tine, set 

A = flow area, ft? 

x - axial distance, ft 

m p = incremental mass flow rate, lb/ft-sec 

u = gas velocity, ft/sec 

g » gravitational constant, ft/sec* 

C p = gas heat capacity, ft/° B 
Tp m flame temperature, °P. 

These equations are applied at each computational increment, as defined by Hie 
mesh for the grain geometry calculation. 

The cross-sectional areas are supplied by the grain geometry calculations. 
The local mass flow Increments, m^, are determined by the grain geometry and 
the local burning rate (sr» Section 5.1 for description of burning rate calculation). 

At each time increment a solution is obtained In the following manner. 

The he-' .-end pressure is estimated, and equations 3-16 to 3-19 are Integrated 
down the port, step by step. When the end of the grain Is reached the nozzle 
choking constraint Is tested. The procedure Is Iterated, with successive head- 
end pressure estimates until the choking constraint is satisfied. The choking 
constraint is applied as follows. First, a throat stagnation pressure Is defined 
as: 

? T * '^KPwT (3 ~ 20) 


r T 

W(HN) 

C* 

AP(NN) 


throat stagnation pressure, lb/ln 2 
mass flow through the throat, lb/sec 
flow characteristic velocity, ft/sec 
throat area, in? 
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Starting with ?„ and a Mach number equal to one at the throat, the lasntropic re- 
lations are used to calculate what the static pressure at the nozzle entrance must 
be 1£ the throat Us to be choked. This pressure, P(NN), Is compared to the static 
pressure resulting from the Integration down the port, P(N). Convergence (choking) 
is realized when these two pressure are within a specified tolerance, the C* used 
In equation 3-20 is the theoretical C* multiplied by an empirical C* efficiency 
(Section 4.1) and an insulation dilution factor (Section S.3). 

The thermochemical properties required for the ballistics calculation ere 
obtained horn the theoretical Performance Module, the nozzle erosion rate re- 
quired to define the instantaneous throat area may be input as a table, or calcu- 
lated internally based on the simple correlations presented in Section S.2. 

3. 3. 2. 3 Other Mass Contributions 

Mass may enter the rocket motor from a gas generator, an igniter, or in- 
sulation oblation. These are treated as lumped, rather than distributed, additions; 
hence, they influence the pressure magnitude but not the shape of the pressure 
distribution down the grain port. For lumped additions, it does not matter whether 
the mass enters at the head end or the aft end of the motor. In constructing the 
Iterations, tt was entered at the aft end (which best represents insulation in general), 
therefore. 


W(NN) = f m e dx + WI + WIGN (3-21) 

or 

W(NN) = W(N) + n Aj fj + WIGN (3-22) 


W(N) = 
Bl 

*1 “ 
r I 

WIGN = 


propellant flow rate 
insulation density 
Insulation exposed area 
insulation ablation rate 
igniter or gas generator flow 
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These masses are Input in tabular form. IgrJter {or gas generated flow la input 
as a function of lime, Insulation density and burning rate are input constants, 
while exposed in* ul at too area is input as a function cf time. Ail of these masses 
are assessed to have tire same thenaochemicaj effects. All are optional end not 
required input. It should not be interred from the Inclusion of an igniter mass flow 
that the analysis can predict ignition delay: it does not do so. The analysis 
ascuaes that burning has started to a motor initially at ambient pressure, and will 
property fill the chamber {also empty the chamber at burnout). Therefore, WIGN 
is an artifice and except where ft Is a significant contribution to the total Impulse 
cf the motor, can be included or excluded as desired. If the user Is seriously 
interested to predicting the ignition transient, the method of Reference (13) 
should be employed. 


3.3 ,2.4 Flow Across Radial Slots (Segmented Motor) 


Pressure drop across a radial slot separating the segments of a segmented 
motor is treated for two cases: (1) where the gas velocity in the slot exceeds the 
gas velocity in the adjacent upstream port: U) where it is less than that velocity. 
In either case, the state properties to the slot are assumed equal to the static 
values cf the gas at the end of the upstream port; slot velocity is then determined 
from continuity in the slot: 


u 


s 


w 


ADD 


(3-23) 


« s = slot gas velocity, to. /sec 

W ADD “ A 0 * rate or propellant from the slot faces, lb/sec 
A s = slot width flow area, in? 
p =• gas density, lb/to? 

It is implied that pressure Is constant in tho slot. Further, it is asserted that 
burning rates are equal over the entire slot face and at the strand rate value. 
Propellant to the slot does not see radiating walls, and to practice the separation 
will be large enough to preclude erosive burning. 

If the entering gas velocity from die slot exceeds the approaching gas 
velocity in the port, the slot pressure drop Is calculated as one Increment with the 
momentum equation. It is as though the slot were physically absent, being filled 
with propellant of high m a » Thus *>i ra' lal momentum Is dissipated, and a 
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psssinu re cbop equivalent to accelerating the added flow to port velocity Is de- 
termined. 

If the entering gee velocity from the slot is less than the approaching 
gas velocity In the port, an empirical relation Is used to find the port area of 
the expanded stream tube over the slot: 

A , = A + (D*/A ) (1 + -f> (3-24) 

pi pm x pm u 


v 

*pm 


stream tube area 

local port area 

slot separation distance 


Knowing the area and static pressure In the stxeamtube, other parameters are de- 
termined by Iteration. The pressure drop across this streamtube Into the down- 
stream port Is determined with the assumption of an Inlet efficiency: 


r l+l 


[(pu*), - { 6 u») l+1 | (3-2S) 


P 1+ j = downstream static pressure 
Pj = upstream static pressure 

ETAIN = Inlet recovery efficiency, typically 0.9 for shaip-edged entrance 
geometry 


Optional capability Is provided to refine the treatment to simulate a slot 
face which Is not perpendicular to the motor axis. Causes of a non-normal face 
Include variable Ignition delays In the slot, grain slump, or tapered Inhibitors. 
The corrections are used to locate the axial position of the outer slot diameter, 
knowing the position of the Inner slot diameter at the grain port. Inputs relating 
to these effects are tables of the amount of burning slot diameter versus time, 
and lateral burn-back versus web. 


The analysis of the slot ballistics Is coupled to the analysis of the slot 
geometry, and overrides the normal grain geometry calculations. A card input keys 
the existence of a slot. 


3. 3. 2.5 Delivered Performance 

The various quantities which characterize delivered performance are also 
calculated with the Grain Design and Ballistics Module. The performance calcu- 
lations within this module do not replace the ones previously described (Sections 
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2.3 and 3.1). They are supplementary; to be referred to only when the Grain 
Design and Ballistics Module Is used in a stand-alone mode. 

Delivered thrust (vacuum) Is given by 

V 

F d = W(N) -Z- „ F + P e A e (3-26) 



W(N) = propellant flow rate, Ib/sec 

V e * exit velocity, ft/sec 

fjp » input nozzle efficiency 

P e = exit static pressure, lb/in? 

= nozzle exit area, in? 

Exit velocity and pressure are determined from the standard Uentmplc re- 
lations. The nozzle efficiency may be determined from the seml-emplrical corre- 
lations, which are amenable to hand calculation. The exact calculations (TD2P, 
etc.) would not be available in the stand-alone mode. As used here, nozzle effi- 
ciency Is an input value. Thrust delivered to the Input, ambient pressure, P & , Is 
given by 


Ambient pressure may be a constant or a table to simulate a flight schedule. 

Total impulse Is calculated as the time integral of Fp . Specific Impulse 

a 

Is calculated on tn Instantaneous basis by dividing Fp by W(H), and on a total 
basis by dividing the Integral of Fp by the Integral of W(N). Insulation weight 
Is not .ncluded because of the convention to base performance on propellant weight 
only; however, both Instantaneous and Integrated Insulation weight flow ere output 
If such correction is desired. 



The loss of performance tn solid rocket motors due to finite rate chemical 
kinetics has been estimated to range between 0.2% for high pressure, low metal 
loading systems, to as high as 5%, or mere, Ic: low pressure, high metal loading 
systems, Bence, if the performance of solid propellant rocket motors Is to be 
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accurately predicted over a wide range of configurations and operating conditions, 
the loss due to chemical kinetics must be accurately modeled. 

The loss of performance which we attribute to finite rate chemical kinetics 
Is usually associated with the difference between the maximum possible sensible 
heat re,ease (chemical equilibrium) and the actual heat release due to chemical 
reactions which occur In the rocket nozzle. The loss Is small If the reactions 
proceed at a rate comparable to a small fraction of the stay time In nozzle. How- 
ever, If the recombination reaction rates are slow all of the heat released by com- 
bustion of the propellant cannot be converted Into kinetic energy. As the ratio of 
nozzle stay time to a characteristic reaction time approaches zero, the so-called 
frozen performance limit Is approached. 

The selection of a method to compute the kinetics loss was not as straight- 
forward as some of the other selection procedures. As pointed out !n Reference 1, 
there are two basic techniques which have been used for calculating finite rate 
performance losses. The first method seeks numerical solutions to the exact 
one dimensional, two phase, reacting gas, equations of motion. The second 
technique is approximate, and Is based upon application of a “sudden freeze" 
criteria. 

Equilibrium flow is assumed up a given point In the nozzle, at which the 
flow Is then frozen. The sudden freeze approximation Is more applicable to liquid 
rockets than solid rocket motors. Such methods were not considered for the pre- 
sent program as their accuracy was not deemed to be high enough. The only ex- 

( 19 ) 

tstlng progr.-n of the first type Is that of Kllegel, et al .As discussed below, 
this program was also considered Inappropriate for the present purpose. 

The one dimensional, two phase, reacting gas program of Reference 19 
yields accurate one dimensional performance predictions, but does not directly 
yield an estimate for the kinetics loss alone. The two phase flow loss Is presently 
calculated on a two dimensional basis using the TD2P program (Section 3.5). 

Thus, If the one dimensional, two phase, reacting gas program Is to be used, the 
two phase and kinetics losses, which are coupled in this analysis, must be separ- 
ated in order to avoid a redundancy. 

An accurate method for separating the two losses could be developed by 
modifying the one dimensional equilibrium program (ODE) to handle finite particle 
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lags. Such an approach was discarded, however, due to the extent of the modi- 
fications required. Other, more approximate means of obtaining a kinetics loss 
could be developed, but not without compromising the resultant accuracy. The 
program described In Reference 19 Is also complex, very large, long running, 
and not easy to use. For all of these reasons alternate means were considered. 

Since the present effort was not to Include any new computer program de- 
velopment, the JANN/iF standard computer program for calculating one dimensional 
kinetics losses In liquid rockets, ODK, was considered. The ODK program Is 
widely used and Is computationally efficient. Complete documentation of the analysis 
and computer program are contained In Reference 5. In Its existing form ODK 
was not applicable, however, since It contained no provision for solid or liquid 
phases. 

The existence of a second phase In the flow can result In significant al- 
teration of the axial distributions of the various flow quantities Cl, B , u, etc.), 
and, hence, impacts the reaction rates and Kinetic losses. The existence of 
finite lags between the particles and gas modifies the effect of two phase flow 
on the kinetics loss, but In a second order manner. 

Various modifications could be made to the ODK program to Include the 
effect of two phase flow on the kinetics loss. In Its unmodified form the ODK 
program numerically Integrates the exact one dimensional, finite rate kinetics, 
equations of motion for a gas only system. A modification of the program to 
allow It to treat the particles <n an exact manner would have resulted In a program 
similar -o that of Reference 19, and would have engendered the same problems. 

A modification based on the adoption of the constant fractional lag approximation 
was also considered. In this approximation the velocity and temperature lags 
of the particles are assumed to be finite, but constant. With this assumption the 
one dimensional, two phase, equations of motion can be reduced to those for an 
"equivalent" gas only system, thereby eliminating the need for solving a separate 
set of particle equations. The required modifications to the ODK program, and 
the ODE program (in order to generate a reference equilibrium solution), were de- 
termined to be substantially less than those required to Implement a solution of 
tho “exact" two phase flow equations. The effort required was, however, still 
greater than desired. 
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Another possible modification, actually a degenerate form of the previous 
one, was also evaluated. In the limit of vanishingly small lags, the constant 
fractional lag equations reduce to the so called equilibrium limit, l.e. particles 
and gas have the same velocity and temperature. The system of equations again 
reduces to that for an "equivalent" gas only flow, but the program modifications 
required to Implement a no lag model wen simpler and less extensive than those 
for the other alternatives examined. While this approach seemed attractive, a 
series of calculations was performed In order to estimate the effect of Ignoring finite 
lags before It was adopted. As a result of these calculations It was determined 
that die zero lag assumption Is quite adequate. It Is esttmeted that the use 
of die zero lag approximation should impact the overall performance prediction 
by less than 0.1* (when the kinetics loss itself Is but a fraction of a percent 
the Impact should be closer to 0.01%). 

In view of the previous discussion, the ODK program modified to treat par- 
ticles on an equilibrium, no lag, basis, was adopted as the method to be used in 
computing the kinetics loss. A kinetic efficiency factor, njau' ls deflned by 
taking die ratio of the I sp (vacuum) calculated by ODK to an equilibrium I gp ob- 
tained by exercising a special option which was incorporated in the ODE program. 
When approached In this manner the kinetics loss is obtained directly; without 
the need to separate out redundant two phase flow losses. The only requirement 
ls that the particles be treated Identically In ODK and ODE, 

Both ODE and the modified ODK program assume the particles and gas are 
In thermal and dynamic equilibrium. The ODE program In Its original form, allows 
for the gas phase to condense as It expands out the nozzle, l.e. the weight 
fraction of the condensed phase may Increase. The modified ODK program, how- 
ever, does not allow mass transfer between phases. : is, to be consistent, and 
to preserve accuracy, a special option was added to the ODE program which allows 
It to "freeze" the concentration of the condensed phases at their chamber values. 
The same dictates of consistency and accuracy governed the treatment of particle 
solidification. The ODE program allows the liquid phase to solidify, and properly 
accounts for the resultant energy release. An analogous treatment of particle so- 
lidification was, therefore, Incorporated Into the ODK program. 

The previously mentioned definition of the kinetic performance efficiency, 
as calculated, may be written In the following symbolic form. 



l 


l 

i 


i 


i 


3-24 






"It! 


“1 


i -i 


1 


akin 


sp ODK 


(3-28) 


»P, 


RE 



I 

i 



wherein 1 aenotes the specific Impulse calculated using the so-called restricted 

sp R£ 

equilibrium option to the ODE program. An additional advantage of using this refer- 
ence method of calculating if. Is that only thermodynamically important species 
sp 

need be considered in the ODK and ODE-RE calculations. Thus, trace species, 
and reactions Involving them, need not be Included in the calculation. 

The following is a brief summary of the ODK program (prior to modification). 

Further details regarding the details of the analysis and the program Itself may be 
found In Volume II of tills report and In Reference * The ODK program calculates 
the inviscid, one- dimensional, nonequilibrium, frozen, and equilibrium (the ODE 
program is incorporated Into ODK), nozzle expansion of gaseous combustion pro- 
ducts. The program can treat large, complex, chemical systems. A maximum 
of ISO distinct chemical reactions and 40 Individual species can be included In 
a given calculation. Chemical reactions are input in stanaard reaction form and 
are translated by an Input processor Into a mathematical form suitable for compu- 
tation. Chemical reactions can be added or deleted by adding (or deleting) a 
single input card per reaction. Up to 10 reactants and 10 products can be Included 
In each reaction. Thera Is also a convenient method for specifying third body 
efficiencies. The coupled set of nonlinear, ordinary differential equations. Is 
integrated using a very efficient, stable, second order. Implicit method. Under 
certain conditions, when efficient to do so, the program switches to an explicit 
Integration method. The program Is designed for englneeiing use. Is well docu- 
mented and user oriented. 

The modifications to the ODK program accomplished during the present 
effort are described In the next section and in Volume II of this report. The reaction 
set and reaction rates, for aluminized propellants, that were employed In the cal- 
culations performed to date are described In Appendix A. The manner tn which 
these reactions and rates were selected Is also discussed therein. i 

3.4.2 Modifications TO The ODK Program I 

! 

Modifications to the ODK program were required tn order to Incorporate 
into It the ability to treat two phase —nansions, with solidification. As discussed 
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In the previous section, particle-gas equilibrium (zero lag) was assumed. Thus, 
the gas and particle velocities and temperatures are taken to be equal. Other 
assumptions which were employed are: 

• There Is no phase change from the gaseous phase to either the liquid or 
solid state. 

•The liquid particles of a given specie solidify at the melting point of the 
solid phase. 

•The particles occupy zero volume and exert no pressure on the gas. 

As would be expected for a program like ODK, a significant effort was re- 
quired in order to identify, program and check out the requisite modifications. As 
mentioned previously, the zero lag assumption allows the two phase flow to be 
treated as an "equivalent" gas only flow. The existing equations employed In 
the ODK program can be made to correspond to the equations for an equilibrium 
(no lag) partlcle-gcs mixture simply by treating the condensed phases „s gas phase 
species with Infinite molecular weights. Thus, the only the equation for the gas 
constant required special treatment. The gas constant Is defined by 


where 


R = 2 CjRj 


R, - — for gas phase species 


( 3 - 29 ) 


R 


1 


R 


= 0 for condensed phase species 


A method for handling particle solidification also had to be Incorporated 
Into the program. Under the zero lag assumption, particle solidification must 
be handled in a special manner, since there Is no explicit heat transfer term to 
govern the rate at which a particle solidifies. It is assumed that when the gas 
reaches the particle solidification temperature die gas temperature remains con- 
stant, l.e., d T/dx =» o, until the latent heat of solidification Is given up by the 
liquid particles. This constraint yields an equation for the rats of change of the 
condensed phase mass fraction, as follows. By noting the constancy of the total 
enthalpy, and applying the restriction dl/dx=0. It Is easily shown that the follow- 
ing equation must be satisfied: 


dT 
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( 3 - 30 ) 
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where N g Is equal 
relations. 


to the number of gas phase species, 

dc s dc z 

"ST = " Hx 


By using the following 

(3-31) 


h -h = aH (heat of fusion) , v 3— 32) 

is rn 


equation (3-30) may be solved for the rate of change of the liquid phase mass 
fraction, yielding, 


dc l 
"3 x 


\ A dC t 

1 „ du . i_ 

■'^m[ U ^ 5 l 1 **\ 


(3-33) 


The ODK program integrates the quasl-ono dimensional equations of motion 
either on a pressure defined, or area defined basis. Therefore, before equation 
(3-33) could be implemented the velocity derivative had to be replaced by either 
a pressure or area derivative. Using the equations of motion given in Reference 
S, it can be shown that equation (3-33) can be written in the following alternate 
forms, for pressure, or area defined problems, respectively. 


~k - j’r 1 T -K - B j pressure defined (3-34) 

-S? = <^npr- | <7# * A, ' B | Area deflned {3 * 3S) 


With dc^dx and dc^/dx defined, the remaining equations can be Integrated. 

This procedure is then applied at each succeeding Integration step until the concen- 
tration of the liquid phase becomes equal to zero (l.e. all of the heat of fusion 
contained in the original mass fraction of liquid has been given up to the gas phase). 

Since this part of the analysis Is new, and cannot be found elsewhere, the 
partial derivatives of the rate of change of the liquid mass fraction, dc^dx (de- 
noted c ), viith respect to all of the other dependent varlalbes are given here, for 
complet<5hes:t. These partial derivatives are required by the Implicit integration 
scheme used In ODK. 
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Pressure Defined Equations 
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Area Defined Equations 
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(3-43) 


where T x = dT/dx = gas temperature derivative In the absence of solidlilcatlon. 

the 0DK computer code has been dimensioned to handle up to 10 condensed 
phase species (or 5 doublets of solid and liquid phase) with 10 discrete melting 
points. This allotment should be adequate to handle any metalized propellant 
system of Interest. 
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While strictly a modification to ODE, and not DDK, the manner In which 
ODE was modified to restrict the condensed phase mass fract.on from changing 
Is described here, as this “restricted equilibrium" modification serves only 
to provide a reference I value for the kinetics loss computation. 

For purposes of computing a reference (or base) I for the kinetics loss 

sp 

calculation an option was Incorporated Into the ODE program which allows it to 
"freeze" the concentrations of the condensed phase species at their chamber values 
(except that corresponding liquid and solid phases are allowed to exchange, l.e., 
solidification Is allowed). 

In order to Implement this option, the ODE program searches the list of 
species to see If aluminum (Al), beryllium (Be), boron (B), lithium (LI) , magnesium 
(Mg), zinc (Zn), zirconium *Zr) , or any compounds containing these elements, are 
present In either the liquid or solid phases. It any such species are found, the 
corresponding liquid and solid pnases aie "enamed, but assigned their usual 
thermodynamic properties. Tats is done at, that the program will not recognize 
them fi'wft-t they are, tnereby restricting these liquid-solid pairs from reacting 
with any species other than themselves. As far as the gaseous species are concer- 
ned, these condensed phases no longer exist, hence, uddltonal amounts of conden- 
sed phase cannot be formed as the flow expands out the nozzle. 

3,5 Two Dimensional Two Phase Module 

3.5.1 '-nproach 

Almost all rocket nozzles of practical interest have a Joss tr, perfonnance 
due to a nor.axlal component of velocity at the nozzle exit plane. This loss is 
usually referred to as a divergence, or two dimensional flow, loss. I f the pro- 
pellant Is metalized, there is an additional loss when the particles cannot main- 
tain a state of dynamic and thermal equilibrium with the gas flow. This loss Is 
referred to as the two phase flow loss. 

The two dimensional and two phase tlow effects are coupled. From the 
discussion In Section 2.4 It Is clour :hat any analytical methodology for predic- 
ting two phase flow loss should be Integrated with the calculation of dl’Xigence 
lass. This requires a two-dimensional flow Held program with two phase capab.llty. 
Both a survey conducted by JAN' . «-e study reported In Reference 1 reco-n- 
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mended, as the best available computer program of this type, the Two Dimensional j 

Two Phase Flew computer program developed by Kllegel and Nickerson^. This 
program was, therefore, adopted for use In the present performance prediction j 

methodology. The basic assumptions employed In the Tti2P program are: \ 

• The flow Is at the choked steady state condition for an axially symmetric ; 

converging - diverging nozzle. 

SThe gas and particle phases each have constant thermodynamic properties. 

•There are no mass or energy losses from the system. i 

«Ihc gas is Inviscld except for Its Interactions with the condensed particles. ] 

•The volume occupied by the condensed particles Is negligible. j 

•The thermal (Brownian) motion of the condensed particles Is negligible. ^ 

•The condensed particles do not Interact. ; 

•the condensed particle size distribution may be approximated by groups 
of different size spheres. i 

•The internal temperature of the condensed particles Is uniform. 

•Energy exchange between the gas and the condensed particles occurs only - 

by convection. 

•The only forces on the condensed particles are viscous drag forces. 

•There Is no mass transfer from the gas to the condensed phase during 
the nozzle expansion. j 

6 

The overall method of solution consists of essentially two distinct parts: i 

a transonic sc’utlon; and a supersonic solution. The transonic solution is required 
in order to provide an Initial line for the supersonic solution. 

The transonic solution I s obtained by an approximate technique. The 
nozzle geometry Is assumed to consist of a conical inlet section Joined smoothly 
to a constant radius of curvature threat. The flow in the conical Inlet portion of 
the nozzle Is assumed lo be a one-dimensional sink flow. The flow in this region 
is coupled to a two dimensional, Sauer type, expansion solution in the throat re- 
gion. The output from this analysis Is a set of flow properties along a supersonic 
start line. The flow field from this initial line out to the nozzle exit is obtained 
numerically using the method of characteristics to sorve the exact two phase flow 
equations of motion. 

This overall procediue works quite well In most cases, however, the naiur® 
of the transonic solution Introduces certain restrictions on the problems which 
can be solved. As previously mentioned the transonic analysis requires that the 
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convergent portion of the mzzle be conical and that the throat be circular, with 
a tingle radius of curvature. Nozzles not of this shape must be “fitted" to conform 
to this geometry in order to obtain a solution. While this can often be accomplished 
without Introducing significant error, the "fitting" procedure Is subjective, which 
is undesirable, and the size of the error Introduced Is net known a priori. The 
transonic method of solution will also fail when the nondinenslonal throat radius 
of curvature (radius of curvature divided by throat radius) Is less than about l.S. 
iind when the conical Inlet angle exceeds approximately 45 s . This fact results in 
two additional restrictions on the n ozzle geometries which may be considered. 

It Is clear from this discussion that a mare exact, less restrictive, method 
for solving the two-dimensional . two-phase flow equations would be desirable. 

With presently available numerical techniques such a method of solution could, 
and should, be developed. 

There is one additional limitation of the TD2P program that has yet to be 
mentioned. If the exhaust nozzle Is contoured In such a manner that the particles 
impinge upon It before reaching the exit plane, the supersonic solution will fail 
terminate) at the first point of Impingement. Xt is a relatively simple matter to 
modify the program to continue the calculation In an artificial manner after Impinge- 
ment occurs. It is not quite so simple, but not too difficult, to Incorporate a 
simple, but reasonable, physical model of Impingement into the program. 

In order to use the TD2P program It Is necessary to determine appropriate 
average values for the perfect gas and particle phase properties. A method for 
obtaining such properties was developed by XI leg el and Nlcherson and was for- 
malized In the Average Gas Properties (ACT) computer program This program 
uses the one -dimensional equilibrium (zero lag) approximation to find a set of 
"equivalent" perfect gas properties such that the solidification area ratio Is Iden- 
tical to the value computed by the ODE program. There are other details involved 
In this computation which will not be discussed at this point. An updated, modi- 
fied, version of the ACT program has been Incorporated In the TD2P program, and 
Is fully described in Volume □ of this report (see subroutine ACT). 

Because the two-phase flow loss Is an Important source cf performance 
loss and Is a function of particle size, the ability to predict particle size will 
have an important Impact upon the success of the computer program over a broad 
range of conditions. The present method for predicting particle size Is described 
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in detail in Section 5.4. The correlation described therein yields the mean particle 
size. The distribution function employed, and the method for selecting the dia- 
meter at each of the particle groups is described in the subroutine DIST write up 
in Volume E. 


■OS 


■.] 


3.5.2 Modifications tr> the TD2P Progr a m 

Numerous modifications and additions to the TD2P program were accomplish- 
ed during the present effort: The more significant of these changes are discussed 
below. 

The TD2P program was originally written In FORTRAN U, and the common 
linkages were completely unsatisfactory for use as a module element. The pro- 
gram had only blank common with multiply defined entries. These were eliminated 
and data storage replaced by multiple labeled common statements. These modi- 
fications were necessary to reduce program size to the 55X requirement for the com- 
bined program. The TD2P over lay structure was revised accordingly. 

The number of options for describing the nozzle wall contour was increased. 

In particular, the ability to use a cubic spline fit contour, rather than a point by 
point specification, was added. A cone, defined by the nozzle end point only, 
option was also incorporated into the program. The maximum number of wall points 
allowed was also increased from 60 to 150 to allow better definition of the contour. 

The TD2F input was revised to be more convenient for the user. Certain 
input items were built into the program, while others were modified or eliminated. 

The partible size correlation presented In Section 5.4 was incorporated 
into the program, as was a subroutine CDISfl to provide a log-normal distribution 
into the specified number of particle groups. 

(23) 

The AGP program was modified and then Incorporated into the TD2P program. 
The AGP program provide!; a rational means for selecting perfect gas properties 
to use in the TD2P calculation. Previously, this p rog ra m had to be executed 
separately and the requisite parameters input to the TD2P program. The AGP pro- 
gram was modified to calculate I on a cne-dlmenslonal, zero lag, basis, to 

sp 

serve as a reference value for the two-dimensional, two-phase loss calculation. 

The manner in which the AGP program calculates specific heats, and the way it 
treats solidification, was also modified. 
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The TD2F prog ram was also modified to provide certain quantities to the { 

roundary Layer Lons Moduli. These quantities are the internally tabulated wall 
geometry, the Mach number distribution along the wall streamline, the two-diaen- ij 

siosa., two-pnase, mass flow rate, and the total pressure and total temperature ' 

at the wall initial 11ns point. The latter two quantities are provided to allow the ]j 

TBL calculated wall pressure and velocity distributions to closely approximate the 
corresponding TD2P value: . 

The object of the HJ2P Module Is, of course, to calculate a performance 
loss efficiency. The manner In which the program war modified to do this Is ] 

described In the next sub-section. \ 

| 

3.5.3 Performance Loss Calculation 

The TD2P Module provides throe types of quantities to the Master Control ‘ 

Module for use in the overall delivered performance calculation. These quantities 
are the nozzle discharge coefficient, C D : and two-dimensional, two-phase, specific 

Impulse, I , and specific Impulse efficiency, nlrwo. as functions of area 

ratio. PTD2P ^ D2P 

The discharge coefficient, C^, is obtained as a direct result of the TD2P ’ 

transonic solution. It is used In obtaining the combustion efficiency from the c* J 

efficiency - equation (3-4), and in computing the ambient pressure correction to 

delivered I - equation (2-4). £ 

sp v 

The specific impulse calculated by the TD2P program I , Is used in 

sp TD2P 

the performance calculation, but never on a absolute basis. At every wall point 

in the TD2P characteristics solution a two-dimensional, two-phase flow, I 

efficiency, tlyQjpif) 1* calculated by ratiolng the TD2P I to a cme-dtoenslonal. 

zero lag, 1 computed at the same area ratio. This one-dlmenslonal, reference 

I jp , Is obtained from the subroutine used by the modified AGP program (AGP was | 

Incorporated into the TD2P program) using the same Ideal gas properties as employed ) 

In die TD2P program. By using a ratio calculated In this manner. Instead of using 5 

the absolute I Jp number from TD2P, the effect of neglecting real gas effects Is j 

minimized. Also, extra cere was taken to accurately calculate the reference I 1 

sp J 

during solidification so as to avoid the use of different assumptions in the two | 

calculations. The equations used in calculating the one dimensional reference I j 

sp 3 

are given below. ] 
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The following equations were iterated using the method of secants to calcu- 
late the ca» dimen si o n al reference I for a given expansion ratio, , . The three 
sen of equations reflect the regions of interest. These are before solidification 
of the condensed phase begins, during solidification, and after solidification ends. 
In the first and last regions, temperature was used as the Independent variable, 
while during solidification, the ratio of the amount of solidified condensed phase 
to total condensed phase was used. 


Before solidification begins; t *. 


where. 
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where. 


a «„ 


“p “g RT 
m 

gas velocity at beginning of solidification 

expansion ratio at beginning of solidification 
weight flow of solid condensed phase 
weight Cow of the gas 

weight flow of the total condensed phase Clqutd plus solid) 
melt temperature of the solid phase 
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The efficiency, *»tx> 2P' 6191x341165 onl y the combined two-dimensional two- 
phase flow effect on performance, as the TD2? and reference I calculations are 
Identical in all other respects. When the nozzle throat erodes, nozzle exit area ratio, 
and, hence, tjjpjP' ls tlm6 dependent. The program calculates a time averaged 
2-D, 2 phase, efficiency In the manner shown In equation (3-3). 


The specific impulse calculated by the TD2P program (as a function of 
area ratio) is also used to compute the erosion loss. While this loss could have 
been estimated based on one-dimensional results (using the ODE program). It was 
felt that a two-dimensional calculation would be superior. As defined herein, the 
erosion loss Is equal to the ratio of the 2-D, 2-phase I at the time averaged 

sp 

expansion ratio, to the specific Impulse calculated at the initial expansion ratio. 


i.e. i 


I 

^erosion = 1 


sp TD2P 

sp TD2P 
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It should be noted that, here again, I values are used on a relative, 

rather than an absolute, basis. This effi8?£icy ls not reported separately, but Is 
combined with the time averaged two-dimensional two phase efficiency, »7 td 2P' 
to yield a single value, ' S6e equMioo (3-2). 

In the present analysis, the fact that the two phase flow loss tends to 
decrease with Increased throat diameter (erosion) was Intentionally Ignored, since 
It was felt that the added loss due to throat roughness would compensate for 
this effect. 


3.6 Boundary Laver Loss Module 
3.6.1 Approach 

IVro dimensional two phase flow solutions Indicate that the bulk of the 
particulate matter gets channeled towards the centerline of the motor due to initial 
effects In the throat region of the nozzle. As a result of this phenomena gas 
phase only boundary layer calculations can usually be made without fear of Intro- 
ducing significant error due to the neglect of the particulate matter. 

The usual boundary layer loss mechanisms which must be accounted for In 


gpww 


rocket engine perfonnance predictions are the energy loss through the. walls due 
to heat transfer, and the momentum loss due to sheer stress at the wall. As dis- 
cussed In Reference 24, these losses can either be found by properly integrating 
the stress tensor along the wall (pressure plus skin friction) or by calculating the 
momentum deficit In the boundary layer at the nozzle exit plane. The latter method 
Is the simplest and easiest to use, and Is die recommended JANNAF procedure 
Following this procedure die boundary layer loss may be calculated as 


i F tt -( 2 »r p U»e cos a ) g 


1 - 


j£_P 

® 0U» 


(3~4S) 


where r, p, U and P are the nozzle radius, density, velocity and pressure at the 
wall exit station, respectively, 6 * and 8 are die boundary layer displacement and 
momentum thicknesses, and a Is the angle between die slope of the nozzle contour 
and the centerline. 


It should be pointed out that this equation gives the complete boundary 
layer loss. The reduction In fluid momentum due to both shear and heat transfer 
losses to the wall is properly accounted for. 

The boundary layer edge conditions that ero required In this formula are 
supplied by the TD2P Program. This allows the effect of two dimensional flow In 
the nozzle to be Incorporated Into the boundary layer solution. This Is the same 
procedure as adopted by JANNAF In the liquid propellant perfonnance methodology 
(l.e. the TDK - TBL Interface). 


Ihe boundary layer thicknesses 5 * and 8 may be calculated In several ways, 
categorized below In order of Increasing accuracy: 


Simple correlational techniques 
Integral methods 

Finite difference and related numerical methods 


The correlational techniques are normally based on semi-emplrlcal cor- 
relations of experimental data for flat plate-zero pressure gradient-flows. In 
many cases the relations are based on Incompressible results and are applied to 
compressible, variable property flows, through die use of compressibility trans- 
formations and reference temperatures (enthalpies). Normally, one of several 
forms of Reynolds analogy Is ap-"<!<5 *o relate skin friction to heat transfer. As 
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a result of their simplicity, these methods have proven to be of great utility in 
preliminary design calculations; however, final engine design calculations are 
normally based on more accurate techniques. 

the majority of rocket engine heat transfer and boundary layer loss cal- 
culations are currently carried out using integral methods such as the TBL program 
described In Reference 7. The integral methods are based upon solutions of various 
integral forms of the differential equations of motion. Ihese methods give two 
ordinary differential equations for die momentum thickness, g, and energy thick- 
ness, j*, (one equation for g only, when the energy equation Is replaced by the 
Crocco relation relating the temperature profiles to velocity, as Is the case In 
TBL). These equations cannot be solved without additional assumptions for the 
form or shape factor, the surface shear stress, r , and heat transfer, q w< 
and the boundary conditions such as wall temperature and pressure gradients. 

The accuracy of Integral methods Is limited by the types of relations 
which are adopted for computing T ^, q w , etc. The majorityof the integral methods 
that have been applied to rocket nozzle flows proceed on the basis of all cr most 
of the following assumptions: 

Tw can be calculated from correlations based on zero pressure gradient 

flows. 

q w can be calculated using Reynolds analogy. 

The shape factor, H=«*/e, can be determined using power law velocity pro- 
files (usually a l/7th power Is taken). 

The Crocco relation can be used to calculate the temperature profiles. 

Reference temperature or enthalpy methods such as Eckert's can be used 

to account for variable properties in the boundary layer. 

All of these assumptions are violated to various degrees in variable pro- 
perty flows In cool wall rocket motors; however. Integral methods of the above 
type. Including TBL, give results which In many cases compare fairly well with 
experimental data. In the supersonic portion of the nozzle, TBL heat transfer 

predictions are generally within 2056 o: the measured values, while in the general 

( 26 ) 

region of the throat, discrepancies of 50% or more have been noticed' The 
large discrepancies between the Integral theories of TBL type and the convergent 
section-throat data can be attributed to the presence of large prossure gradients 
In the throat region which seriously loiote some of the assumptions upon which 



1 
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the analysis Is based, the occurence of "lamlnarlsatton* of the boundary layer 
due to die pressure gradients, and die sensitivity of the calculations tz. assump- 
tions made about the boundary layer flow In the chamber. 


While of major Importance to local heat transfer calculations these dls- j 

crepancles in die throat region are not serious from the standpoint of performance 
calculations. For most rocket motors T8L, or related Integral methods, should ’ 

be able to calculate boundary layer losses to within about 20%. Considering 
that boundary layer losses are usually about 1%, and rarely exceed 2% (small 
motors, of short firing duration being one of the exceptions), the total impact of j 

such a 20% error on overall performance loss should be restricted to several tenths .5 

of a percent; an amount that should prove tolerable In most cases. 

There are more accurate methods avail able for calculating boundary layer 

performance loss, such as the Mass Addition & 'indary Layer Program^) and the 
(28) 

BUMP Program' , These methods, are based on finite difference or other num- 
erical solutions of the turbulent boundary layer partial differential equations. 

While the advances achieved In recent years allow such solutions to be achieved 

In reasonable amounts of computer time, the trade-off between lnceresed accuracy 

and Increased operational complexity (as well as run time) did not appear to favor 

the uae of a finite difference type approach In the SPP Program. This conclusion • 

follows from the usually restricted size of the boundary layer loss, which logically 

precludes a msjor effort to calculate It. 

The previous discussion outlines why It is felt that Integral methods, and 
TBL In o-tlcular, are adequate for the task at hand. More accurate methods such 
as MABL and BUMP are not warranted for general solid rocket motor performance 
prediction. However, in certain cases the Increased generality of programs like 
these may bo the only way -o obtain an adequate prediction for boundary layer per- 
formance loss. If such cases can be Identified, a program like MABL can serve, 
not as an Integral pert of the performance prediction methodology, but as an adden- 
dum to It; to be used as required. On the other hand, there will be classes of motors 
for which data and or correlations exist which will adequately be able to characterize 
the boundary layer loss. This will be especially true for the class of motors for 
which boundary layer losses are small. In such cases the use of any analytical 
boundary layer solution may be obviated. 
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The SPP Program has been structured to account for all of these possibilities. 
The TBL program was selected as the standard method for computing the boundary 
layer loss, and was Incorporated into the program as a module. Two other methods 
for detenu inlng the boundary layer loss have been explicitly provided for. A simplified 
semi-empirical boundary layer loss correlation can, on option, be utilized in lieu 
of the TBL program. This correlation is described in Section 4.2.4, The boundary 
layer loss can also be directly input to the SPP Program. In such a case, the input 
value of th8 loss is used in deference to all others. In this manner, should the 
circumstances warrant it, a more exact boundary layer method may bo utilized to 
calculate the boundary loss in an external, uncoupled, fashion. The computed 
value may then be Input to the SPP program. 

The TBL program does not allow for a direct assessment of the time de- 
pendent nature of the boundary layer loss. Por motors of all but very short firing 
duration, this is of little consequence to die accuracy of the overall performance 
prediction. Since wall temperature Is an input to the TEL program, the variation 
of the boundary layer loss with wall temperature (and therefore with time) may bo 
evaluated by performing 5 series of boundary layer solutions wherein the nozzle wall 
temperature distribution is varied. To perform this calculation in a rational man- 
ner a thermal analyzer type, transient, heat conduction cods should be externally 
coupled to the boundary layer solutions to provide the required wall temperature 
distributions. 

3.6.2 Modifications to the TBL Program 

Very ilttla modification was required to Implement the TBL program as a 
module In the SPP program. Most of the modifications were related to the Incor- 
poration of additional labeled commons to allow TEL to communicate with other 
modules, as appropriate. Many of the Items that must normally be Input to TBL 
are Instead obtained through this linkage. For most problems, the need to specify 
the bulk of the remaining Input items was obviated by presett.og them in program 
data statements . Thus, very little user specified Input Is normally required to 
execute TBL as an element of the SPP performance prediction methodology. 

The only other modification of substance to TBL was the incorporation of a 
calculation of the boundary layer induced performance loss decrement, ^1 
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This quantity was calculated as follows: 






sp. 


TBL 


A F BL 
m TD2P 
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where jF^ Is given by equation (3-4S), and ls the no2zle mass flow rate 

calculated In the TD2P module. 

The TBL calculation Is carried out at the average chamber pressure, and 
the initial expansion ratio. The boundary layer loss is assumed to be Independent 
of erosion Induced effects for the present purposes. 


3.7 Inter-Module Interfaces 

In order for the various programs contained within the SPP program to function 
together, automatically, the program modules must be able to communicate with 
each other, as required. A perfunctory description of the Inter-module Interfaces 
ls given here to provide the reader with an overall Idea of how the elements of the 
SPP program are linked together. More detailed descriptions of the program linkage 
structure may be found In Volumes II and HI of this report. 

The following five simplified diagrams are designed to Indicate only the 
nature of the Information which flows out of the five basic computational modules, 
and Its ultimate destination. In actuality, all of the Inter-module communications 
are processed through the Master Control Module. Each block In these diagrams 
represents a module. The following abbreviations are used for convenience: 

Master Control Module: MC 
Theoretical Performance Module: ODE 
Grain Design and Ballistics Module: GD&BC 
Xlnetlcs Module: ODK 
Two-Dimensional Two-Phase Module: TD2P 
Boundary Layer Module: TBL 
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4. PERFORMANCE LOSS CORRELATIONS 



4.1 Combustion Efficiency 
4.1-1 Approach 

The combustion efficiency of aluminized solid propellants Is ordinarily close 
to 100 percent. Departures from Ideal characteristic velocity generally result 
from wall losses (as can exist with an exposed steel chamber In subscale test 
hardware), from particle lags In the throat region, and from dilution by ablatives 
(generally corrected out In data reduction). These thr> . , 'fects will be further 
commented upon subsequently. Exceptions to the ef '<*> . combustion of aluminum 
are found with low burning rate propellants, operation at low chamber pressure, 
and tests In combustors of exceptionally low L*. Interest In spacecraft motors, 
upper-stage motors and sustalner operation has provided substantial data showing 
the loss of combustion efficiency at low chamber pressures and burning rates^ 4 ^. 

Therefore, the correlation approach began with a theory that alum', am combus- 
tion efficiency Is a function of propellant burning rate. Burning rate bears directly 
on those processes which Influence the competition between extensive metal 
agglomeration and rapid metal Ignition at the propellant burning surface^ 1 ^ . Burning 
rate directly influences the combustion zone temperature gradient, the particle 
heating rate, i— 1 the gas velocity normal to the propellant surface. Reported 
effects of propellant lormuiatlon ore relatable to changes In burning rate^ 4, ^ 
cs long a' ‘he oxidizer and binder provide enough energy for me.al Ignition (always 
true in ale ..i.ilzed propellants ci oractlca! Interest, not always true for other metals). 
Likewise, reported effects of pressure are relatable to burning ‘ate. The motor L* 
does not appeal to be Important so long as It Is greater than approximately 50 In., 
according to experiments measuring complete combustion dlrectly(16),the Implication 
Is that complete combustion is largely determined by what happens at or near the 
propellant surface. Motor C* efficiency data were accumulated Ir several groups 
for the purpose of constructing, verifying and refining the correlation with burning 
rate, and finally to demonstrate the completed correlation on an aprlorl basis. 

Statistical analysis Is a usef .l diagnostic tool to provide Insight Into the 
nature of poorly-understood phenomena, but has proven misleading when used to 
generate correlations for refined predictions. The present computr-r program devel- 
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opmeot was predicated upon the premise that factors affecting the specific impulse 
efficiency of aluminized propellants are not poorly understood, and that existing 
theories and models can be evaluate* and combined tc achieve the desired accuracy. 
Further, statistical analysis is not warranted when the data base is insufficient 
compared to the potential scope of application. While die amount of BAXES motor 
data Is considerable, tee fact teat it is a standard test motor becomes restrictive; 
potential governing parameters will vary together in a corresponding manner and so 
cannot be meaningfully separated. It ts also true teat tee great majority of the 
BATES data are at a standard test condition, and cover a fairly narrow range of bum 
ing rates and aluminum concentrations. Therefore, tee use of statistical analysis 
was limited to a supplementary role. 

Data groups that were studied are listed in Table 1. 

4.1.2 The Basic Correlation with Burning Kate 

The basic correlation of C* efficiency with burning rate is shown In Figure 4-1. 
These data consist of Groups I and II, excluding data for propellants having aluminum 
concentrations below 10 percent or greater than 24 percent. Each point represents 
an average for a series of flrii,gs of a given propellant, in a given motor, at a par- 
ticular test condition. Where BATES data for a given propellant existed in g re ps, 
each group was averaged separately. It is observed that th bulk of tho data exist 
between burning rates of 0.3-0. 6 in./sec. and, within teat range, the correlation 
Is not meaningful at first glance. It is only when a broad rar ,*e of burning rates 
are examined that tee correlation appears significant. 

It is observed teat C* efficiency falls rapidly below burning rates of 0.25 
in./sec. This effect persists In tee full-scale motors which are induced in tnat 
range. Therefore, C* efficiency Is not expected to depend upon scale. The pro- 
blem with the steep behavior at low burning rate Is teat the uncertainty In predict- 
ability becomes large In that regime. At high burning rate, the efficiency appears 
to approach 100 percent. 

4.1.3 BATES Motor Reproducibility 

The basic correlation with burning rate was examined in tee light of tee repro- 
ducibility of the BATES motor C* efficiency data, particularly in view of the above- 
mentioned confinement of most of the data. 


i 
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Table 4-1 



DATA USED TO DEVELOP, REFINE AND VERIFY 
C* EFFICIENCY CORRELATION 


GROUP I DATA - CONSTRUCTION OF CORRELATION 


a) Propellants In 6 x 11.4 Rohm & Baas Motors 

LPC-343-1401 

LPC-343-1402 

LPC-343-1403 

LPC-58QA 

LPC-S&2A 

LPC-586A 

I.T-C-62SC 

NOTS X-32 

RHP-112 

TPH-1011 

LPC-36S-1SAL (NAS 1-9182) 

NASA/Langley (NAS1-10956) 


b) Full-Scale Motors 


ASTROBEE-F (Susialner) 

Aerojet 


AVANH 

Lockheed 

Company records 

ATS (Etvi- Sumer) 

NASA/JPL 

Ref ( 29 ) 

ATS (Internal- Burner) 

nasa/jpl 

Ref ( 29 ) 

HARPOON (Lockheed Testwelght) 

Lockheed 

Company records 

1S6-S 

Lockheed 

Ref ( 30 ) 


GROUP H DATA - BATES MOTORS, VERIFICATION AND REFINEMENT 1 
a) Conventional Propellants 


Ref ( 17 ) 

Ref (17) 

Ref ( 17 ) 
Company Records 


Ref (14) 
Ref ( 18) 


RHP-112 

LPC-S43A 

RHP- 113 

LPC-543D 

RHP-161 

LPC-580A 

RHP- 162 

RDS-S01 

RHP-163 

RDS-507 

TFG-3016D 

RDS-510 

TPH-10RS 

UTF-11475 

TPH-1011 

CSSA 

TPH-1066 

VBB 

TPH-8163 

SPIS-31 

AHP-2969 

SPIS-32 

ANB-306F 


ANB-3105 


AAB-331B 




Table 4-1 (Continued) 


b) Advanced or Experimental Propellant* 


ARCHT-S015 

RHU-105 


ARCHT-9018 

DLH-258 


ARCHT-9021 

ANB-3406 


ARCHT-9024 

ANP-3413 


ARCHT-9027 

ARCEF-147 


ARCHT-9030 

LXA-100 


UTT-1S908 

SPIS-26 


UTX-14617 

SP1S-27 


UTX-14682 

UTX-14631 

UTP-1S151 

VLU 

VLZ-IJ 

VMO 

' m DATA - DEMONSTRATION OF CORRELATION 


2.75 FFAR 

Aerojet 

Ref (31) 

2.75 FFAR 

Hercules 

Ref (32) 

SRAM 

Lockheed 

Ref (33) 

HYDAC 

Lockheed 

Ref (34) 

HI-EX (Experimental) 

Lockheed 

Ref ( 35) 

MINUTEMAN III 3rd 

Aerojet 

Ref (36) 

FW-4S 

UTC 

Ref (37) 

SURVEYOR RETRO 

Thlokol 

Ref i 38) 

ANTARES I 

Hercules 

Ref ( 39 ) 

ANTARES H 

Hercules 

Ref ( 40 ) 

SPARROW 

Rocketdyne 

Ref (41) 

SAM-D 

Thlokol 

Ref ( 42 ) 

-°ARTAN-2nd 

Thlokol 

Ref ( 43 ) 

156-5 

Lockheed 

Ref (30) 


FIGURE 4-1 BASIC CORRELATION OF C EFFICIENCY 



BURNING RATE (IN'SEC) 




Within the range of burning rates 0. 3-0.6 In/sec, the correlation exhibits 
a maximum deviation of 1.2 percent and an average deviation of 0.6 percent. All 
groups of SAKS data were examined for reproducibility. The maximum deviation 
within a group was found to be as high as 2.0 percent, and the average deviation 
within a group was found to be as high as 1 .4 percent. These results were not 
flash-ln-the-pan. One group consisting of eight nominally Identical firings exhibited 
an average deviation of 1.2 percent, and average deviations covering all groups 
were typically greater than 0.5 percent. The accumulated average deviation was 
0.8 percent. Thus, the accuracy of the correlation could be said to be comparable 
to the reproducibility of the data. This was not particularly satisfying because 
unsound correlations might turn out just as good within this range of burning rates, 
forcing reliance on the data outside of this range to support the proposed correlation. 

Some attention was devoted to the reason for the data scatter In the BATES 
motors. The principal uncertainty In the C* measurement Is caused by nozzle erosion; 
C* requires relating the measured pressure Integral with an average throat area, 
a factor not present In the specific Impulse measurement. All other factors. In- 
cluding Instrumentation and propellant weight, are also present In the specific 
Impulse measurement. Average deviations In specific impulse efficiency ranged 
from 0. 1-0.3 percent, clearly and significantly superior to the C* error. Reported 
nozzle erosion rates, and throat averaging, are based upon pre-fire and post-fire 
measurements, they do not reveal the nature of the throat behavior in the course 
of the firing. A cause of variability In throat erosion behavior Is the repetitive use 
of the same nozzle Insert for economy. It was confirmed by AFRPL that this was the 
general practice and, within a given group of data, one can observe this taking place 
from the Initial throat sizes. The effect Is aggravated In high erosion situations. 
Therefore, unless this data scatter In C* Is acceptable. It would be recommended 
that fresh or noneroding nozzles be used on each test for Improved accuracy. 

The Improved accuracy of the specific Impulse measurement afforded an 
alternative opportunity to look for burning rate effects through specific Impulse. 

4.1.4 Effect of Burning Rate on Specific Impulse Efficiency 

Additional confirmation of the effect of burning rate on performance appears 
In the specific Impulse efficiency data, which are more accurate than the C* efficiency 
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data, and, therefore, lend more credence to the effect. This Is shown for Group 
□a propellants In Figure 4-2. Specific Impulse efficiency Is plotted versus aluminum 
concentration, as Is often done. Each datum point here Is an average of an over 
all series of tests for each propellant. The number placed next to the datum point 
Is the average burning rate. Note that, at each aluminum concentration. Increasing 
burning rate Increases specific Impulse efficiency. This Is believed to reflect C* 
efficiency changes. Therefore, wide variations in burning rate can explain why 
specific Impulse efficiency cannot cleanly be plotted versus aluminum concentration. 

The darkened points In Figure 4-2 denote one type of propellant In which 
aluminum concentration was systematically changed, with burning rate approximately 
constant. All of the of the other Ingredients are the same. Thus, the darkened 
points would basically reflect the expected change In nozzle efficiency with Increa- 
sing aluminum concentration. 

Assuming the burning rate correlation to be valid guide to C* efficiency, 
the nozzle efficiency was then computed from this correlation and the specific 
impulse efficiency data. Results are plotted versus aluminum concentration In 
Figure 4-3. Note that the spread of data points Is considerably less In Figure 
4-3 than In Figure 4-2. This reflects normalization of the influence of burning 
rate. 

4.1.5 Effect of Binder on Specific Impulse Efficiency 

A curious oddity appears In Figure 4-3. Except for the propellant In which 
aluminum - Mcentratlon was systematically varied. It cannot be concluded that 
efficiency decreases from 15 percent aluminum to 19 percent aluminum. Rather, 
the data Indicate that efficiency Is higher at 16-18 percent aluminum than at 15 
percent. Upon closer examination of these conditions. It was found that the anomaly 
could be explained by a binder effect. Data from 15 percent to 19 percent aluminum 
could be systematically organized by binder, ar.d the more efficient binders were 
concldentally at the higher aluminum concentrations. Within the Group Ila pro- 
pellants, the order In the direction of increasing efficiency was polyurethane, NC/ 
NG, CTP3 and HTPB, PBAA and PBAN. 

This binder effect was Incorporated as a correction on the correlation of C* 
efficiency with burning rate. By correcting C* efficiency for binder, the Figure 
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FIGURE 4-2 CORRELATION OF SPECIFIC IMPULSE EFFICIENCY WITH 



PER-CENT ALUMINUM 



FIGURE 4-3 CORRELATION OF SPECIFIC IMPULSE EFFICIENCY WITH 
ALUMINUM CONCENTRATION, NORMALIZED FOR BURNING RATE 
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4-3 data should collapse Into place and clearly show the expected dependence 
of specific Impulse or nozzle efficiency on aluminum content. This is shown in 
Figure 4-4, which incorporates all Group n data and similar corrections for Group 
nb binders. 

The value of the C* correlation In specific impulse efficiency correlation 
is summarized as follows: 

a) If specific impulse efficiency in BATES motors is correlated i .iy with 
respect to aluminum concentration, the standard deviation Is 0.8 per- 
cent and the probability of accuracy within 2 percent is 94 percent. 

b) If the burning rate correlation is taken Into account, the standard de- 
viation is 0.4 percent and the probability of accuracy within 2 percent 
is better than 99 percent. 

c) If both binder and burning rate effects are taken into account, the stan- 
dard deviation is 0.2 percent, and the probability of accuracy within 

2 percent is better than 99 percent; this Is comparable to the reproduci- 
bility of the specific Impulse data. 

Thus, the BATES motors were successfully correlated within the limits of 
the data received. 

4.1.6 Effects of Extremes of Aluminum Content on C* Efficiency 

The data correlations do not Include propellants having aluminum concentra- 
tions in excess of 24 percent. Such propellants were tested in BATES motors as a 
matter / academic interest, but are too fuel-rich to be of practical interest. As 
might be expected, values of C* efficiency fell below the Figure 4-1 correlation, 
metal agglomeration was undoubtedly extensive. On the other hand, values of 
nozzle efficiency increased. The Increase In nozzle efficiency could be explained 
by a reduced total particle concentration resulting from the fuel-richeness. In any 
event, it was concluded that should propellants, of practical Interest, having high 
aluminum concentration ever be developed, they would probably obey the correlation. 
Therefore, the correlation is not modified for high aluminum content. 

C* efficiency data for low burning rate propellants having aluminum con- 
centrations below 10 percent fell above the Figure 4-1 correlation. Aside from this 
fact, a correction would be warranted by the premise that combustion inefficiency 
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FIGURE 4-4 CORRELATION OF SPECIFIC IMPULSE EFFICIENCY WITH 
ALUMINUM CONCENTRATION, NORMALIZED FOR BURNING RATE AND BINDER 




stems from the aluminum; therefore, some means should be afforded to approach 
the unmetalllzed propellant condition. A correction factor In the form of an equa- 
tion was added to accomplish this. 



% 
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4.1.7 Effect of Scale on C* Efficiency 

The C* correlation should be unaffected by scale Inasmuch as Group I 
data Include large motors. However, there Is one factor related to scale which 
does need to be addressed. A large motor is likely to be fully Insulated, whereas 

a subscale test motor (Including BATES motors) Is likely to have exposed walls. 

(44) 

Measurement of heat loss In BATES motors Indicates that a 1 percent loss may 
be assigned to the chamber. A fully Insulated motor would have negligible heat 
loss because the thermal energy in the insulation is recovered by ablation until 
the last Instant of the firing. Therefore, In constructing the final correlation, all 
subscale data were shifted upwar J by one percent, and a one percent loss factor was 
Included In the correlation for an uninsulated motor. 

4.1.8 Treatment of Discharge Coefficient Effect 

In addition to combustion and heat losses, C* efficiency also embodies 
a two phase two dimensional discharge coefficient effect. The effect of C p on 
performance is calculated within the TD2P module, so the danger that a performance 
loss would be counted twice In the specific impulse calculation was presented. 

To assure that this would not happen, the following approach was taken. 

!he value of C* efficiency determined from the correlation was accepted 
as valid, being based upon motor data. In the overall performance calculation the 
C D related portion of the C* efficiency Is Isolated, and removed. The remainder 
Is a combustion efficiency factor (and heat loss). This procedure is detailed in 
Section 3.1. In the Grain Design and Ballistics Module the C* efficiency Is used 
directly In the calculation of chamber pressure. 

4.1.9 Verification of Co relation 

The final C* efficiency correlation consists of a table look-up, with cor- 
rection factors In equation form, and Is presenter ir> Table 4-2. A comparison with 
motor C* efficiency data from Group III is presented in Table 4-3 . 
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Table 4-2 


C* EFFICIENCY FUNCTION 


”C* = [ K + TT (10 ° ' 10 ] x b x 0 

percent aluminum; a = a for a £ 10, a •- 10 for a > 10 
1.00 for fully Insulated motor, 0,990 for uninsulated motor 
binder constant 


NF, C = 1.003 
PGA/NC, C 1.008 
PBAA, C = 1.006 
PBAN, C = 1.006 
HTPB, C = 1.003 
CTPB, C = 1.000 
NC, C = 0.998 
PU, C = 0.992 


burning rate constant 


k 

r 

k 




SI. 4 

0.50 

98.6 

93.1 

0.60 

98.9 

94.0 

0.70 

99.1 

94.6 

0.80 

99.2 

95.1 

0.90 

99.3 

95.6 

1.00 

99.4 

96.0 

1.20 

99.6 

96.4 

1.40 

99.7 

96.7 

1.60 

99.8 

97.0 

1.80 

99.9 

97.7 

98.2 

a 2.00 

100.0 


'1 

i 


>100.0, = 100.C 


Table 4-3 
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Group III motors were selected bacause of certain features bearing upon 
specific impulse predictability. In the C* context, thev also afford a wide range of 
burning rates, aluminum concentrations and hinder types. The range of variability 
of 3.3 percent lr. C* efficiency provides a reasonable basts for testing the method; 
clearly, the test has little meaning If the range Is comparable to the accuracy. 

It Is observed that the average deviation of the predictions Is 0.4 percent. 
Although the predictions are uniformly distributed between high and low values, 
the overpredlctions tend to be larger In magnitude. This would be consistent with 
the uncertainty of accounting for insulation dilution m the actual data. Insulation 
dilution was not considered in these predictions {this effect Is discussed In Sec- 
tion S.3). Because these motors were not used In developing the correlation, 
and the accuracy Is a small fraction of the range of variability. It Is concluded 
that the Table 4-2 correlation is satisfactory. 

4.2 Specific Impulse Efficiency 
4.2.1 Approach 

The developme it of simple expressions for each source of specific Impulse 
loss affords a means for the rapid and convenient assessment of motor delivered 
performance. Such expressions therefore provide an Important and useful supple- 
ment to the analyr-col computer program subroutines which deal with these losses 
In a more exact and detailed mat er. 

In order to enhance the scope of applicability of the simplified method, a 
seim-anai . ral rather than statistical correlative approach was adootec. The dis- 
advantages of the statistic,.’ approach were discussed In t' c context of combus- 
tion efficiency correlations. In correlating specific impulse efficiency, there is 
the added consideration that tt is a function or n. merous, Interrelated processes. 
Depending upon the propellant and the motor design, different physical mechanisms 
may dominate the loss. Therefore, correlations of delivered specific Impulse or 
total nozzle efficiency In terms of alam.num concentration and scale parameters 
will be deficient to the extent that the .ndlvldu3l loss components do not vary 
systematica’' with these variables. Ac. appirac’' that addresses *he Individual 
loss components is preferred, and a sec i-analytlcal approach Is bettei able to 
d~ so. 


1 i 
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The resulting expressions were derived from analytical or seml-emplrlcal 
expressions appearing in the literature, or by fitting the results of the analytical 
computer programs used In the SPF program. A group of motors was selected that 
would serve to highlight various loss components, and furnish a basis for their 
verification. The final expressions were compared to the measured total loss for 
all of these motors . 

The motors selected for this exercise, and the reasons for their selection, 
are listed In Table 4-4. These motors represent extremes of the motor design vari- 
ables that influence delivered performance. The noz zle efficiencies of the produc- 
tion motors are plotted versus aluminum concentration In Figure 4-5, and are com- 
pared to the Figure 4-4 BATES correlation. It .s observed that there is no correla- 
tion with aluminum content. It is also Interesting to note that most of the produc- 
tion motors deliver lower performance than the BATES motor even though most are 
larger. Therefore, factors other than aluminum content, or size per se, must be im- 
portant. These factors do not vary appreciably In BATES motors used as a standard 
propellant reference, but must be considered in predicting motor performance In 
general. The factors which distinguish production motors are contained within the 
expressions for the loss components. 

The seml-emplrlcal efficiencies described In Sections 4,2.2 through 4.2.7 

are In terms or per cent nozzle thrust coefficient, C p , losses Instead of fractional 

I efficiencies, 

sp 

4.2.2 Divergence Loss 

The simple half -angle correction, modified for contoured nozzles^ * , is 
employ Thus, the details of the source flow from the throat region, and the 
coupling of the two-dimensional particle lags, are not taken Into accojnt here. 

This Is a potential source cf error in the simplified analysis. There Is insufficient 
Information to isolate these coupled effects, either experimentally, or from para- 
metric TD2P solutions, to create a Justifiable simplified expression. One dlt .ulty 
In performing a parametric analyses Is that the full shape of the nozzle has an Impor- 
tant Influence on the coupling of the two-dimensional and two-phase flow effects. 

The expression used for divergence loss is as follows: 

or + 4 tv 1 

r DFv - * r ’0 1 - cos { ^>i (4-1) 

3 = half angle of conic i. nozzle. Included angle of contoured nozzle 

Spy 3 exit angle of co-.i . i -ozzle 
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FIGURE 4-5 ABSENCE OF PRODUCTION MOTOR PERFORMANCE CORRELATION 












4.2.3 Kinetics Loss 

The kinetics loss expression originally proposed simply stated that this 
loss is 1/3 of the difference between theoretical shifting and theoretical frozen 
performance at the particular expansion ratio. It came from a study of various pro- 
pulsion systems, including liquid engines and fluorlnated propellants under condi- 

(45) 

tions that would emphasU^ this loss . However, subsequent comparison with 
results of the kinetics calculations from the ODK program indicated that a pressure 
correction should be provided to reduce the predicted loss at higher pressures. 

The final form of the kinetics loss expression is therefore: 


, = 33.3 1 - 


theoretical frozen Isp 
theoretical shifting Isp 


,200 , (4-2) 

• '*7~ T>200psl 


A thenr.ochemical calculation is necessary to exercise this expression, and care 
should be taken that the frozen result is for the same expansion ratio as the shifting 
result (some computer programs furnish optimum results at different expansion 
ratios). The difference between shifting and frozen performance will vary with 
propellant chemistry, pressure and expansion ratio. 


4.2.4 Boundary Layer Loss 

Boundary layer loss is expressed in the form of a time-dependent heat transfer 
expression, with a correction term for nozzle expansion ratio: 


’Ibi =c i'^t | 1 * 2exp (-c 2 p°-Vd°- 2 )J[ 


1 + 0.016 (( - 9) j (4-3) 


P = pressure, psi 
D ( = throat diameter, in 
t * time, sec 
( = expansion ratio 

Ordinary Kozzle : 


Steel Nozzle: 


Cj = 0.00365 
C 2 = 0.0009J7 

C. * 0.00506 

c 2 * 0 
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The time-dependence Is In exponential form to represent transient heat-up. 
Wall heating with time Is an Important factor to consider In motors of short bum 
time {e.g., less than 4 seconds). The heat transfer variables are combined in the 
well-known form for pipe flow. The dependence on expansion ratio represents 
the effect of increased nozzle surface area. 

The bracket term containing the time-dependence states that the heat trans- 
fer with a cold nozzle (t=0) will be a factor of 3 greater than with a fully heated 
nozzle (t -«*>). This factor comes from comparing a typical equilibrium wall temp- 
erature-gas temperature difference with the difference between ambient temperature 
and the gas temperature. It is consistent with the differences between hot wall 
and cold wall heat losses derived from computer calculations, and plotted in Ref. 

1 . The time constant (C^) comes from an analysis of the transient heating of 

a standard BATES motor. The coefficient (Cj) was Judiciously selected based upon 
a study of the following: a direct measurement of heat loss in a BATES motor (Ref. 

44 ); computer calculations of boundary layer loss (Ref. 1 ); measurement and 
analysis of nozzle heat transfer coefficients (Refs. 46,47); deductions of the 
boundary layer losses In standard 3KS and 1GKS motors by knowing the total measured 
loss with unmetalllzea propellant (Ref. 48,49) and assuming that divergence loss and 
kinetics loss are correctly described by the equations above. 

The constant preceding the expansion ratio term was deduced from the HI-EX 
motors as lollows. The motor data revealed a decrease m nozzle efficiency with 
increased expansion ratio. A series of TD2P calculations provided results showing 
an Increase In nozzle efficiency (decrease In 2- D lag losses) at high and Increasing 
expan. . , n ratios*. Assuming the expression for kinetics loss to be correct, the 
remainder >■ as assigned to differences In boundary layer loss with expansion ratio. 

A second set of constants is provided for a "steel nozzle," which Is primarily 
addressed to the Hercules 2.75 motor. These constants assume a constant cold 
wall (Cj-O) , and a higher temperature potential (higher Cj). They should be 
limited to a motor employing a solid steel nozzle of relatively thick walls. 

*An Important practical application of this result Is that propellants of higher 
metal content and total solids ma, .< considered for high expansion ratio than 
the optimum at low or standard e>-ans!on ratios. 
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4.2,5 Ttoo-Phase How loss 

The expression for two-phase flow loss Is more complicated, attempting to 
take Into account different regimes of concentration effects and Interdependencies 
of particle size and throat size effects. The form of the expression and the values 
of the various constants were obtained by fitting detailed solutions appearing in 
References 1, 50-52, and from parametric TD2P solutions obtained during the pre- 
sent effort. 

The two-phase flow loss Is given by: 


”TP a ^.1S ( 0.06 D Cg 

£ = mole fraction of condensed phase, moles/100 gm 

D p = partie’e size, u 
P = pressure, psl 
r = expansion ratio 
D = throat diameter. In. 


If r*0.09, C 4 = 

O.S 










D t '• 1: C 3 = 

9.0, 

C 5 

= 

1.0, 

C 6 

= 

1.0 




1 < D { <2: C 3 = 

9.0, 

C 5 


1.0, 

C- 

o 

= 

0.8 




D t -.2: and Dp < 4: 


C 3 

= 

13.4, 

c s 

= 

0.8, 


= 

0.8 

4 s Dp <8: 


C 3 

= 

10.2, 

C 5 

= 

0,8, 

t 


0.4 

D p .P: 


=3 

= 

7.58, 

S 

= 

O.S. 

C 6 


0.33 

If {<0.09, = 

1.0 










D t <l: C 3 = 

30.0, 

~5 

- 

1.0, 

C 6 

= 

1.0 




l<D t <2: C 3 = 

30.0, 

C 5 

- 

1.0, 

C, 

o 

S 

0.8 




Dj •, 2 and Dp < 4: 


C 3 


44.5, 

C 5 

= 

0.8, 

C 5 

= 

0.8 

4<D p -s8: 


S 

= 

34.0, 

S 

e 

0.8, 

C 6 


0.4 

Dp -8: 


C 3 

= 

25.2, 

C 5 

= 

0.8, 

C 6 

- 

0.33 


D p = 0.454 P 1/3 r 1/3 1 1 - exp 
L* - motor characteristic length. In, 


-0.004L*) 1+0 


,0451^ j 
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Hie expression branches first as a function of exhaust condensibles. It 
then branches as a function of throat diameter, and then with particle diameter, for 
large throat diameters. In general, the dependencies upon particle size and throat 
size Increase as these quantities get smaller. There Is a slight dependence upon 
expansion ratio. Effects of pressure, exhaust condensibles, throat size and particle 
size are coupled. In that particle size Is a function of the other three; the result Is 
that the loss will Increase slightly with pressure, strongly with exhaust condensibles, 
and go through some orttmum respecting throat size. 

The expression tor particle size combines theories of particle growth by 
condensation in the chamber and collisions in the nozzle, as discussed In Sec. 5.4. 

4.2.6 Submergence boss 

The expression for submergence loss Is In the form suggested by Kordtg and 
Fuller J , with constants modified based upon the data appearing In Ref. ( 1 ). 

P, 0.8 c0.4 

hSUB^ 0 - 0584 ZfJ (4 ' 6) 

“t 

P « pressure, pst 

£ = mole fraction of condensed phase, moles/100 gra 

A* = nozzle entrance area/nozzle throat area 
S = length of subT.crgcr.ee/iength cf internal motor 
D t = throat diameter, in- 

4.2.7 bosses Not Corslderec 

Expressions were not Included for the following: flow separation In highly 
over -expanded nozzles; shock losses stemming from rough nozzle surfaces caused 
by nor, uniform erosion oi ^palling, expansion ratio losses stemming from erosion or 
nonopttmum expansion by design, "the last Item Is not really a loss, and may be 
accommodated by comparing predicted and measured performance at the actual 
overage expansion ratio. It has s.gntficgace lit spriorl prediction where the actual 
throat erosion u not Known. Ar-el/sos oi shock losses and nonuniform erosion are 
not yet veli-dsveiopce tor curler, rrtdlctron. Separa.iec requires considerable over- 
expansion, not uxoiy except dur , portions of Ignition and tall-off which are 
generally :econpeqvonttai awe : ’ - total Imculse. 






4.2.8 Verification of Simplified Loss Predictions 


J 


Results of applying the simplified methodology to the selected gropu of 
motors are listed In Table 4-5. The simplified performance predictions were ob- 
tained In the following manner. The Individual nozzle loss per cent efficiencies 
given by equations (4-1) to (4-4), and (4 S) were combined to yield an overall 
nozzle Cp efficiency, 

^CF 100 - (^bl + "DIV + n f3N r n SUB + 

In performing this calculation average values for pressure, expansion ratio and 
L* had to be estimated. Particle size was obtained from equation (4-5). A C* 
efficiency was obtained from Table 4-2, and combined with the efficiency 
to yield an overall delivered I efficiency, 

T)j = t? c * tj cf (4-8) 

sp 

The measured values of I efficiency were obtained by dividing the measured 
sp 

delivered specific Impulse (delivered to ambient pressure) by the theoretical 
specific impulse (at ambient back pressure) calculated at the estimated average 
expansion ratio. 

It is observed that all predictions are within l.S percent of data, and that 
the average deviation Is 0.5 percent over a range of variability of 5.9 percent. 

The overpredlctions and underpredictions of efficiency appear to be equally distri- 
buted In number, but the overpredlctlans tend to be greater in magnitude. 

Th SATES motor results are added to the production motor results In order 
to comment upon the performance of NF propellants relative to standard propellants. 
These results indicate that the supeilor performance of the NF propellants In 3ATES 
motors was an artifice of the particular propellants and the BATES motor. The NF 
propellants were of substantially higher burning rate. Improving C* efficiency Ir. 
addition to the Inherent Improvement from the binder itself. The higher burning rate 
necessitated a larger throat size, which In Itself unproved the ncczle efflclen y 
1 percent due to lower two-phase flow losses. Thus, the Indtcatea 2 percent 
Improvement would not be prea'lctec to be maintained m a full-scale motor where 
bum rate and throat size would be fixed by design. The Inherent Improvement 
In combustion efficiency could also ue compromised by an increased kinetics 
loss with fluorine at high expansion rv.:-. 
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Range of variability = 5.9% 
Average deviation = 0.5% 


S. OTHER CONSIDERATIONS 
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5.1 Burning Rate 

5.1.1 Approach 

Motor burning rates cannot, as yet, be predicted on a strictly a priori 
basis. Therefore, burning rate predictions must be based on experimentally mea- 
sured strand burning rates for the propellant of interest. The problem of calculating 
motor burning rate then reduces to accounting for the factors which cause It to 
differ from the strand buring rate 

At a given set of conditions (i.e. pressure, temperature) strand burning 
rate Is a constant. Motor burning rapes are design dependent, however, being a 
function of the local convective and radiative heating at the burning siiface^' 

The differences between motor burning rates and strand homing rates are often large 
enough to Impact pressure and thrust to a greater extent than the desired tolerance 
level of five percent; hence, these differences must be accounted for. The sophis- 
tication of the analysis used to convert strand data to motor bum rate Is dictated 
by a compromise between accuracy and the information readily available for an 
a priori calculation. 

One cf the more comprehensive coupled analyses of motor burning rate 
appears in Sef.rence 13. That analysts was constructed for the special pijposes 
of predicting Ignition delay. Ignition spikes, the success or failure of command 
termination, and delivered Impulse during transient start-up and shut-down lr, 
multiple su,k,-restatt operation. For those purposes, a comprehensive combustion 
and energy exchange analysis was required, Including an In depth heat transfer 
solution for the solid phase. The accurate prediction of steady-state ballistics 
followed as a by-product of the accurate prediction of transient ballistics. The 
disadvantages of this analysis are Its complexity, size and run time. The run 
time and size, with a limited number of axial Increments, are comparable to 
goals established for the entire Grain Design and Ballistics Module; coupling It to 
the much finer grid required by the gram design solution was considered imprac- 
tical for a priori predictions. Therefore, the use of this analysis should be 
limited to the detailed study cf transient ballistics for which It was created. 

The present program Is chiefly ccncor *»d v ith steady or quasi-steady ballistics. 
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The approach that was taken was to construct a closed-iorm expression 
consisting of an erosive term, a radiation term, and a transient term: 


r = r Q (1 + fir) (I +r‘) 


r = strand rate. In. /sec. 
o 

fir = steady-state corrections 
r' = transient correction 


5.1.2 Erosive Burning 

An extensive review of erosive burning was performed by LPC In associa- 
tion with the development of a computer program to predict the ballistics of nozzle- 
19 54) 

less rocket motors ' .All of the erosive burning laws were found to be deficient 
In some respect, but the best and most comprehensive equation appears to be that 
of Lenoir and Roblllard: 

c G 0 ’ 8 e 0 r 

ar - e *p <- - V — > < s - 2 > 

Cj, c 2 = constants 

G = local port mass flux, lb/ln? -sec. 

x = local axial position. In. 

o = propellant density, lb/ln“ 


In this v ork, the r appearing In the exponential w as changed to r Q to remove 
the iterative step: It Is a minor correction, within the uncertainty In the values of 
Cj anc. The dimensionless c 2 was left at 53 per Lenolr-Roblllard. The constant 
Cj may be expressed In analytical form, but then includes other corstants of more 
or less uncertainty. The value of this constant was theiefore established from 
selected motor data wherein erosive burning would be expected to be significant: 
0.0068 In Inch"’ 8 - pound"’ 8 - second’ 8 units. 


5.1.3 Radiation Contributions 

Radiation can be a significant source of combustion augmentation in large 
(13) 

motors . In the case of unmet -tlilzed propellant, radiation will be imparted by 
hot nozzle and Insulation surface; as well as by the hot gases flowing inside 
the motor. In the case of metal,.. «.l propellant, the particle-gas cloud shields 
the propellant from the hot inc. >ut is Itself a source of significant radia- 

tion heating . 
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Radiation contributions to motor burning rate have been studied using the 
sophisticated transient analysis computer programs of reference (13) which dealt with 
motors containing unalumlnl^ed propellant and Reference (55) dealt with motors 
containing aluminized propellant. Closed form expressions are constructed as 
ratios of radiation flux to combustion flux: a) in the case of unaluminized propellant. 




r F 

0.1 

1 ^ 

1.125 T* f - 

)(Aj/A b ) + 

(A^/Ag) 

+ 

l-exp(-xP) 

(Tp - T* ) 


(5-3) 


b) in the case of aluminized propellant 

[l-exp(-\ P}] (Tp - Tj) 
= * 1-8 O s C~r-J F 


(5-4) 


where: 


= Stefan- Boltzmann constant, 3.52 x 10 


-14 


BTU 


ln. a -sec- 6 K* 

average view factor from walls, based on Ref. (13) computations 
thermochemical flame temperature, °K 

propellant surface temperature, °K=850 for AP composite propel- 
lants, 523 for double-base propellants 

fourth power of ratio of insulation wall temperature to propellant 
flame temperature, taken from typical Ref (13) output 

exposed insulation area, in? 

propellant burn surface area, in? 

projected annular nozzle area, radiating surface, in? 


0.1 

T F 
T s 

0.125 

A- 

a n 

0.410 = nozzle analogy to 0. 1 25 above 
= emlsslvlty, (lb. /in?) * 


X 

P 

C s 

1.8 


BTU 


pressure, lb/ln? 
propellant heat capacity, 

conversion factor 


and 


X = 


/14.69 


(5-5) 


0.001 + 0.04 (Percent Alumlnum)/100 
from Reference (55). 

Results obtained for motors / heroin the radiation contribution was signi- 
ficant Indicates that equations (5- *' ' - r 0 yield good results, without having 

to apply any additional "correction - 
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5,1.4 Transient Term 


Effects of P on Instantaneous burning rate were approximated using the ex- 
pression derived by Paul and Von Elbe: 

r ‘ = 27=*P (5-6) 

o 

n = pressure 

a = thermal diffuslvity, ln?/sec 
P = rate of change of pressure 

5.1.5 Other Effects 

Effects of acceleration on burning rate, brought about by spin-stabilization 
of motors, were not included in the program. If desired, such effects can be in- 
corporated at a later date. 

Burning rate augmentation along bondlines has been encountered in certain 
motors. The phenomenon has been attributed to propellant nonuniformity. No 
attempt was made to account for it because It is to be avoided in the future. 

Recent analyses have begun to address two-dimensional erosive burning. 

In these analyses bum rate differences between the tips and valleys of stars and 
keyhole slots are accounted for in addition to longitudinal differences. Although 
progress a' ing these lines is encouraging, the incorporation of two-dimensional 
erosive burning into the program would have required further new developments, 
and, therefore, was beyond the scope of the present effort. 

Conditioning temperatuie is accounted for through program input, l.e., 
the input strand burning rates should correspond to the conditioning temperature 
of interest. 

5.1.6 Strand Data as Input 

In recognition of the fact that the burning rate pressure exponent, n. Is 
not a constant, and that conditioning temperature effects vary with pressure and 
temperature, It was decided to Inn.;: strand burning rate data in tabular form. 
Logarithmic interpolation Is used within this table. 
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Best Available Copy 




S . 2 . 1 Approach 

Changes In nozzle throat area with ttir.e mainly Impact the chamber pressure. 
Thrust and specific Impulse will be affected to a lesser extent, but the effect can 
become significant If the throat area changes appreciably. 

Because nozzle erosion Is a time-dependent aerothermochemlcal effect. It 
Is not sufficient to correlate erosion by steady-state heat transfer, parameters 
which would Include pressure, throat size and temperature. The time-dependence 
of nozzle heating rates must be accounted for, particularly for short bum time, 

A parameter related to the oxidation potential of the combustion products Is also 
required. Combustion of contemporary propellant formulations produces exhaust 
gases which contain several chemical species that are reactive with the free carbon 
contained In throat materials. It follows that these additional considerations will 
vary with the formulation of the propellant and the type of throat material. 

The prediction of erosion rates In carbonaceous and graphitic throat materials 
has been the subject of several computer program developments combining chemical 
kinetic theory and heat transfer. These efforts are very sophisticated and Include 
experimental work as well as theoretical analysis. However, when applied to motor 
firings, accuracies could net confidently be expected to be better than 30 percent 
without a posteriori adjustment of constants. 

According to the original program plan the nozzle erosion rate was to be 
specified by means of a user determined Input table. During the development of 
the computer program It was agreed that provision for some means of calculating 
the erosion rate would be desirable, since not all potential users could be expected 
to have the capability to Independently determine It. Rather than burden the SPP 
program with additional complexity. It was decided to Incorporate an erosion rate 
capability on a seml-emplrlcal basis. Towards this end, the CMA^^ and ASTHMA^ 7 ^ 
computer programs were run for a matrix of conditions consisting of two propellants 
(zero aluminum and 16 percent aluminum), three pressures (500 psi, 1000 psl, 

2000 psl), and three materials (carbon phenolic, ATT graphite and pyrolyt'c graphite). 
For each material the calculated erosion rates were fit as a tabulated function of 
temperature and parameters characterizing heat transfer and combustion products. 
Temperature, in turn, Is derived from t v e propellant flame temperature and e tlme- 
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dependent heat transfer expression. Final adjustment of constants followed 
application to motors representing each material. The exact nature of the erosion 
rate as a function of time has beon retained. Thus, when warranted, the present 
correlations may be overridden by data, or the results of another analysis, deemed 
more appropriate for the problem at hand. 


5,2.2 Form of the Nozzle Erosion Correlation 

Nozzle erosion rate was correlated In the following form: 



^ - i; a c c h 

(5-7) 

where is 

a time dependent film coefficient having the form. 


c,/c • 0.8 „ 

c h - tct hr ) ~rar 

D t * M 

(5-8) 

c 3 

= a constant. Including conversion factors 


D t 

= throat diameter. In. 


m 

= propellant weight flow rate, lb/sec 


\ 

= throat area In? 


X 

= gas thermal conductivity, BTU/ln.-sec.-°R 


4 

= gas viscosity, lb/ln,-sec 



= gas heat capacity, ft/"R 


* t 

= erosion rate, in. /sec 


°rp 

= material density, lb/ft 3 


8 C 

= dimensionless erosion rate parameter 


C H 

= heat transfer parameter, lb/ft s -sec 


12 

= conversion factor 


For each of the three materials considered, equation (5-7) yields the relation. 

\ = 

0.173 ftj, Cjj Carbon Phenolic 

(5-9) 

II 

0.109 S' c C H ATI Qraphlte 

(5-10) 

't = 

0.086 8J. C H P>Tolytlc Graphite 

(5-11) 
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(56,57) 

Hie Aerotherm computer programs calculate g' c as a function of wall 

temperature, for a given propellant and material. The results were found to be 
dependent on temperature and the corrosivity of the propellant combustion products, 
but were essentially Independent of pressure. Thus, the erosion parameter ^ , was 
correlated In the following form, 

fi' c = c 4 (D l (3) (5-12) 

where c^(T) Is a tabular function of temperature, and f(8) is a function of the pro- 
pellant products corrosivity parameter, 9. The temperature in equation (5-12) 
should be wall temperature, however, wall temperature was converted to flame temp- 
erature by assuming the equilibrium wall temperature to be 70 percent of the flame 
temperature. The flame temperature is obtained from a thermochemlcal equilibrium 
solution (ODE program). The time dependence of the wall temperature, expressed 
as a flame temperature. Is obtained from a standard form of transient heather “re- 
pression: 

T g = T o + < T f * T 0 > j 1 * <>XP -C S ( Y ,0 ’ 8t <5 ‘ 13) 

At time equal to zero the temperature is equal to the initial temperature, T o P>K), 
while at longer times, the temperature approaches the propellant flame tempera- 
ture, TpPK) . The constant c^ is found by linearly Interpolating In a stored table, 
with T^ as the argument. There Is a separate table for each of the three materials 
considered. 


Th 

meter 9: 


-orroslvlty of the combustion products is characterized by the para- 


M 5 



Mg = molecular weight of gaseous combustion products 

M q = molecular weight of :..utal oxide 

A - concentration of metal oxide, moles/100 gm 

^ th 

X, = mole fraction concentration of the 1 oxidizing species. 

1 There arc 5 oxidizing species considered, H.O, CO., 
1/2 O z , C, OH Z 2 
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The parameter 8 Is computed by the ODE program and transmitted to the Grain 
Design and Ballistics Module, where It Is used. Based on the matrix of calcu- 
lations that were carried out, the following relations were obtained for f(fl) (see 
eq. 5-12): 

f(g) = 1 Carbon Phenolic (5-15) 

j 2 (no dependence on 0) 

f(8) = (^-) ' ATJ Graphite (5-16) 

f(3) ~ (j ^24 Pyrolytic Graphite (5-17) 


Thus, to summarize, for each material the erosion rate, r t . Is found from the one 
of the equations (5-9) - (5-11). C H Is given by equation (5-8) and 0 c by equation 
(5-12) (with the subsldary relations (5-13) - (5-17)). The constants c^/Cp, In 
equation (5-8), and c $ in equation (5-13) were adjusted based on comparisons with 
a limited amount of erosion data for each of the three materials. The final values 
selected for these constants are: 


c 3 /c p = 1.682 
c $ = 3.368 


(5-18) 


Other nozzle materials can be added to the computer program by repeating the pre- 
viously described procedure for the materials of interest. 

For operating pressures below 300 psl, the phenomenon of nozzle deposition 
rather than erosion Is likely to be encountered In aluminized propellants. Depo- 
sition -z not been studied to the degree of erosion. Limited data suggest a 
rate of -0.001 In. /sec. for every 50 psl below 300 psl. This is offered as a guide 
in constructing an Input table for use with low-pressure motors. 


Insulation Degradation 


5.3.1 Approach 


Insulation has a minor effect on the Interior ballistics of motors. However, 
it Is enough to warrant accosting for this effect, but not enough to conduct de- 
tailed analyses. Therefore, a simple methodology was adopted. 
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S.3.2 Performance Degradation 
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From the standpoint of program accuracy and efficiency multiple thermo- 
chemical calculations, as the Insulation becomes an Increasing fraction of the 
total mass flow, would not be Justified. Instead, one matrix of thermochemical 
calculations was performed varying aluminum content from 0 to 25% and Insula- 
tlon/propellant ratio from 0 to 0.1. The following expressions were fit to the 
results: 

-jjv = ll -[0.5 - 0.025 (%A1)] R (5-19) 

0 = |l for %A1 i 20 

1 -[ 0.7 - 0.035(%A1)] R (5-20) 

1 for %A1 > 20 
T 

= ll - 1.9 R (5-21) 

F o = lo.Sl for Ri 0.1 

where R denotes the lnsulatlon/propellant weight flow ratio, and the subscript o 
refers to values with no Insulation. The performance degradation due to Insulation 
Is considered only In the Grain Design and Palllstlcs Module and is not transmitted 
to the Master Control Module. 

5.3.3 Determination of the Insulation Weight Flow Ratio 

Coupling the Insulation exposure to the grain regression pattern would not 
be Justifies ■>/ the magnitude of the Insulation effect. Rather, a table Is provided 
for a priori Input of exposed area versus time. Ablation rate and Insulation density 
are single Inputs. Therefore, at any time, the flow rate of Insulation Is determined 
and, by comparison with the propellant flow rate, the value of the ratio is established. 
The thermochemical parameters used m the ballistics analysis are then degraded 
accordingly. 

5.4 Aluminum Oxide Particle Size 
5.4.1 Approach 

The two-phase flow loss Is an Important source of performance loss and 
Is a strong function of particle size, * the ability to predict particle size will 
have an Important Impact upon the s„: e. ,1 the computer program over a broad 
range of conditions. 
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Considerable particle size data have been reported In the literature. Al- 
though the collection methods and the relevance of the results have been open to 
question (58), the existing data must be accepted in lieu of anything better. A 
number of mathematical models of oxide particle formation and growth have also 
been published In the literature. 

The approach that was taken was to construct a seml-emplrlcal expression 
for the mean particle size. The form of the expression Is analytical in origin, 
ano the constants result from a fit of data covering a wide range of conditions. 

A log-normal distribution Is constructed about the mean, and selected sizes repre- 
senting distribution Intervals are used in the TD2P module. 

5.4.2 Experimental Data and Pure Empiricism 

(59) 

Aluminum oxide particles were first collected and measured by Brown 
and Sehgal^®®! Brown's work provided the first Indication of their rounh-order-or- 
magnltude size, and Sehgal showed the size to be a function of pressure. Cheung 
and Cohen^ 61 * extended Sehgal's work to show dependendlces no*, only upon pressure 
but also upon residence time and aluminum concentration; no effect of propellant 
aluminum particle size was detected. Since that time, the Air Force has Inves- 
tigated particle size In a variety of motors which has formed a basis for Air 

( 44 ) 

Force performance Investigations . Having acquired data over a wide range of 
motor sizes, the Air Force chose a simple correlation based on nozzle throat diameter. 

The problem with a correlation based on throat diameter is that It does not 
addrer 'he controlling physical phenomena. This is a drawback of empirical 
performance correlations based on size parameters, in general. The reason that 
such correlations work in an approximate sense is that many of the parameters 
tend to vary together (o.g. , throat size, motor 1 * , propellant weight or weight 
flow, bum time, oil tend to Increase in larger motors). However, when Improved 
accuracy Is desired. It is necessary to accommodate variables outside of the cor- 
relation range Including combinations of these variables which are exceptions to 
general trends (e.g., a decrease Ir. L* with increasing motor size) and therefore 
sources of inconsistency or error. Because the scope of computer program appli- 
cation may encompass a wide range of variables, singularly and in combination, 
the throat size correlation would not be expected to produce the desired accuracy. 
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The throat diameter correlation was constructed from data over a relatively 
narrow range of pressure and aluminum content, but over a wide range of sizes, 
wherein L* and throat diameter tended to vary together. Effects of pressure and 
particle concentration, as well as L*, must be Included If the scope of propellants 
and operating conditions Is to be adequately covered. Cheung and Cohen actually 
found a decrease In particle size with Increasing throat diameter In their data 
because they increased throat diameter In order to reduce pressure (and average 
L*) In fixed hardware. Similarly, nozzle efficiencies reported from BATES motor 
data^ 44 '^'®^, again In a fixed size, tend to infer a smaller particle size with 
decreasing pressure (increased throat size and smaller mean L*). The BATES 
data also indicate a smaller particle size with decreasing aluminum content, and 
with changes In binder type that would result In decreasing Al^O^ concentration. 

The phenomena of oxide particle formation and growth must be addressed In some 
manner In order tc account for these detailed effects. 

S.4.3 Models and Semi Empiricism 

Mathematical models of oxide particle growth have been reviewed by Brown ^ S5 ^. 
The Cheung and Cohen paper Included a model describing particle growth by con- 
densation and agglomeration, assuming that all particle growth takes place within 
the motor cavity. They used this model to correlate their data and explain the 
trends observed. They then extrapolated to large boosters and predicted the sizes 
actually measured from these motors two years later. The principal effect was 
residence *t.re, relatable to the motor L*. The model also predicts particle size 
to be dependent upon pressure to the one-third power and Al^Ch, concentration to 
the one-third power. The residence time dependence Is more complicated because 
condensation and agglomeration have different functional time-dependencies. 

Assuming that these time dependencies can be combined In a single term represen- 
ted by L* to a power, there results the following seml-emplrical expression: 


i 


D p ~ P 1 '" 3 $ 1/3 t* m (5-22) 

This expression was fit to the Air Force data, and resulted In a very good 
correlation as shown In Figure 5-1 (m=0.8). The correlation was superior to one 
simply based upon throat diameter. However, when other data were considered. 
Including lag loss Implications In large motors having small nozzles, the L* de- 
pendence was concluded to be too . . 
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Another class of models assumes that particle growth occurs chiefly In the 
nozzle by collisions resulting from velocity gradients. Such a model would not 

have an Independence, nor does It address pressure and concentration effects 

( 66 ) 

from the standpoint of physical chemistry. The model of Crowe and Wllloughly' ' 

Is representative of this class. A seml-emplrlcal correlation of the data was con- 
structed based on this model, and was found to be poor. Therefore, It was con- 
cluded that a model of particle growth should not emphasize growth In the nozzle. 
This Is fortunate from the standpoint of the complexity of the lag loss analysis. 

Finally, a correlative expression was constructed combining the elements 
of the Cheung-Cohen and Crowe -Willoughby analyses. This assumes that particle 
growth occurs in both the chamber and the nozzle, but each to a lesser extent than 
predicted by either theory alone. It also serves the purpose of reducing the In- 
dependence of growth In the chamber under the circumstances where It was found 
to be excessive. The expression fit to the data Is as follows: 

D p = 0.454 P 1/3 e 1/3 [l'- exp(-0.004L*)] (1+0. 045 D { ) (5-23) 

D p In microns 
P in lb/ln? 

£ In moies/100 gm 
L* In Inches 

D ( = throat diameter. In Inches 

Thus, a limiting diameter Is reached In the chamoer at large 1.*, depending upon 
pressure c.‘ i concentration, and additional grow*h occurs In the nozzle depending 
upon throat size. A comparison with the data Is shown tabulated In Table 5-1. 

The average deviation is 18 percent, which represents a variance In specific Impulse 
efficiency of less than 0.5 percent. 

Additional studies of particle formation and growth are in progress In the 
context of particle damping of combustion Instability. Observed trends support 
the theory that substantial growth must occur In the motor cavity. Results of 
these studies may be useful In future modification of the above correlation. 
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Table S-l 

COMPARISON OF PREDICTED AND MEASURED PARTICLE SIZES 



Range of variability = 12. 8 U 
Average deviation = 17.7% 


5-14 










6. COMPARISONS WITH MOTOR FIRING DATA 




While the SPP program was, by design, based almost entirely on pre- 
existing, proven, analytical and semi-empirlcal techniques, there was no 
guarantee that the Integration of all of these methods was properly implemented. 
Nor could It be known, ahead of time, tf the combined techniques would be 
capable of meeting the accuracy goals which had been set (i.e., c-llvered 
Isp to + 2 %, delivered thrust and total impulse to + S%): for a wide range of 
motors , propellants and operating conditions . 

In order to validate the present computer program, a series of compari- 
sons with full scale motor firing data was carried out. To date, four program-data 
comparisons have oeen completed. The first three comparisons were basically 
used as vehicles to check out the operation of the computer program. The three 
check out cases were selected to provide a reasonable range of motor sizes, pro- 
pellants and operating conditions. The motors considered are listed below, 
together with some of the characteristics that led to their selection . 

Thlokol Extended Delta : complicated three dimensional grain, 
relatively large size, high expansion ratio. 

AeiQlet 2.75: small, short bum time, multiple nozzles, high 
throt t erosion. 

Lockheed bRAM : end burner, unusual nozzle shape, high 
pressure, low expansion ratio. 

As a result of their diversity, these motors provided a reasonable 
thorough basis for evaluating the performance of the program. The fourth com- 
parison with motcr data was canied out after the previously mentioned comparisons 
were completed. This was a so-called sealed envelope comparison. The Air 
Force selected the motor to be considered and transmitted to us only the infor- 
mation required to execute the program. The motor chosen for this comparison 
was the C-4 Third Stage. It represents a motor of much current Interest, having 
a high energy propellant, complex grain geometry and reasonably high expansion 
ratio. The completely a priori performance predictions for this motor were then 
delivered to AFRPL and evaluated. 

The results of each of these four comparisons are Individually discussed 
be'ow. The results for all of these cases are then summarized in Section 6.5. 
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6.1 Extended Delta Motor 

The present results for the Extended Delta motor were obtained by 
executing the complete performance prediction methodology on a completely 
coupled basis. For reasons discussed subsequently, this was the only one 
of the first three check-out cases for which this was done. 

The theoretically calculated delivered specific Impulse for this motor 
was found to be 0.46% less than the measured value. The simplified semi- empi- 
rical performance prediction correlations under predicted the measured value by 
0.79%. The results of the ballistics calculations are discussed below. 

The predicted pressure-time curve for the Extended Delta motor is 
compared to the actual pressure trace in Figure 6-1. It is observed that the 
qualitative and quantitative features of the data are well-predicted. Pressure 
is slightly underpredicted for the initial portion, and slightly overpredicted 
for the latter portion. The action time Is somewhat overoredlcted . 

The predicted action time average pressure was 534 psia , compared 
to the actual result of 545 psia, for an error of 2.0%. Average thrust was 
underpri dieted by 3.1%. The action time was overpredicted by 4.3%, in part 
because of a propellant weight error attributed to the fact that the web 
dimerslons had to be measured from a motor drawing. Correcting for the 
pressure and assigned web enrors, burning rate was underpredicted by 2.4%. 

This burning rate error .n part explains the low initial pressure anu longer 
action time. The pressure becomes high near the end of the burn In part 
because the throat erosion was underpredicted by 18.7% and in part because 
of the icreased input web. As a result of compensating effects (i.e., thrust 
underpredlction, weight flow overprediction) the predicted total impulse was 
only 0.5% above the measured value. 

A burning rate error of 2.4%, if uncompensated by errors in the other 
parameters governing motor pressure, would produce a 3.4% pressure error for 
a 0,3 exponent and a 6.1% pressure error for a 0.6 exponent. Therefore, an 
error even in the course of strand burning rate measurement could make it 
difficult to predict pressure to within 5%. Sources of compensating (or aggra- 
vating) error would be in the grain geometry, nozzle erosion and C* efficiency 
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computations. That there Is a grain geometry erroi Is evident lr. Figure 6-1 In that 
the predicted pressure dip precedes the actual pressure dip even though the 
predicted burning rate was low. It follows that the predicted pressure dip 
should have occurred later than actual. However, the prediction of the basic 
ciiaracter of the pressure curve shows that the grain design program Is properly 
computing the geometry In general. In this case, the error In the throat erosion 
partially compensated for the error In the burning rate. 

(13l 

6.2 SRAM First Pulse 


The primary reason for selecting SRAM as a validation case was to 
provide a means to check out the grain geometry and ballistics computations 
for an end-bumer. The motor also has a unique nozzle geometry and r Igh 
operating pressure. A full SRAM firing consists of two pulses. Calculations 
were carried out only for the first pulse. 

The current grain design analysis does not contain a model for *he meniscus 
burning phenomenon encountered In the SRAM motor. Thu"., the ballistics pre- 
dictions for the SRAM motor had to be corrected for the > A meniscus burning. 
Neglecting meniscus burning, the computed pressure and thrust levels were 
expected to be. and were, low. The computer output was corrected for meniscus 
burning by raising the pressure at each instant of time and by reducing the value 
ol each time in accord with the Increase in effective burning rate. Thus, the 
fact that a given K N is encountered sooner, was accounted for. The basis for 
the pressure correction at each instant was the previously known effect of the 
phenomenon on the grain regression. Although the important result was the 
successful operation of the computer program for an end-bumer, the corrected 
pressure-ttrne trace is shown compared with data in Figure 6-2. 

The shape of the pressure trace is well-predicted. The predicted action 
time average pressure was low by 1.5%. The predicted motor burning rate as 
corrected was low by 1.3%, and was compensated by an underprediction of throat 
erosion of 2.5%. Action time average thrust was underpredicted by 0 . 8 %. The 
underprediction of pressure was compensated by a 0.9% overprediction of 
delivered specific Impulse, which is well-within the program goal. 

The specific Impulse prediction for this case was not based on the motor 
operating conditions output by theGraln Design and Ballistics Module, due to the 
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neglect of the previously mentioned me: "is burning phenomenon. Average motor 
pressure and erosion rate, based on the data were Input to the program, together 
with an average motor L*, and the ODE, TD2P and TBL Modules were executed. 

The Kinetics Module was not utilized since the loss is Insignificant at the SRAM 
operating pressure (the simplified kinetics loss correlation predicts a 0.1% kinetics 
loss for this case) . 

6.3 Aerojet 2.75 FFAR( 31) 



The ballistics prediction for the Aerojet 2.75 inch motor (Ref 31) 
encountered an anomaly. As shown in Figure 6-3,, there Is no resemblance 
between the predicted and actual pressure-time trace. A disturbing circumstance 
was an oscillograph trace for a ”2.75 Inch motor" which had been previously 
furnished to Ultrasystems by AFRPL, and which showed progressive burning. 

Tils piece of data will be discussed further subsequently. In view of this 
startling result, an Investigation was conducted as to the nature of the problem. 

After checking the inputs to the computer program, a discussion was held 
with members of the LPC Engineering Department who had 2 . 75 inch motor exper- 
ience. It was reported that there is no reasonable way in which a circular-port 
grain could show a neutral pressure trace, and that inhibiting the aft portions 
would have little effect on the progressivity. 

Nt rtheless. Ref (31) was checked against the engineering drawing of 
the grain and there appeared a basis for suggesting a modification to the Input 
description of the Inhibited surfaces. A test case was run with this modified 
input, but showed little change in the progresslvlty. However, because a 
different ansi ' did result, it was possible to conclude that the grain geometry 
subroutine was a t least iecognizlng inhibited surfaces and probably correctly 
accounting for them. 

The next step was a conversation with Aerojet. It was established that 
the neutral character of the pressure trace was reproducible over an extensive 
test history. The explanation offered was a combination of super erosive burning 
and considerable nozzle throat erosion. However, this explanation proved 
unsatisfactory for the following reasons: 
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Figure 6-3 Predicted and Actual Pressure Traces For Aerojet 2.75-Inch Motor 
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a) The computer program adequately predicted the reported measured 
throat erosion (11.3% error) . It was verified that a single nozzle 
had not mistakenly been utilized In the calculation (the motor has 
four nozzles) . 

b) An approximate calculation established that the rate of throat 
erosion would have to be Increased by a factor of 2 In order 
to produce a neutral pressure trace. Discussion with LPC 
engineers Indicated that a different throat material, such as 
an asbestos-based material, would be required to produce 
such a rate. 

c) The computer program closely predicted the reported measured 
motor burning rate, when corrected for the difference between 
predicted and measured average pressure (1.4% error). 

Therefore, the factor increase between motor and strand rates 
was satisfactorily predicted. 

d) The Initial port-to-throat diameter ratio of the motor is 2 , which 
is too high a value for exceptional erosive burning. Moreover, 
the aft end port, which would be expected to exhibit the highest 
eroslvlty , was Inhibited. 

e) An approximate calculation established that, to produce the 
measured average burning rate and a neutral pressure trace with 
the measured throat erosion, the initial erosive burning would 
have to exceed that produced In nozzleless rocket motors and 
the burning rate G-dependence would have to be equal to the 
pressure-dependence (like a hybrid rocket fuel containing some 
AP). Both of these requirements are unrealistic. 

f) The long tall-off exhibited by the motor Is indicative of some 

anomaly. In seeking to explain the anomaly, discussions were 
held with LPC engineers regarding the possibility of unbonds. 
Inhibitor failure or cracks. An unbond was considered unlikely 
in view of the encapsulation of the propellant in this particular 
design. Inhibitor failure would not significantly alter the pro- 
gresslvlty, as mentioned earlier. The possibility of a reproducible 
cracking of the grain was then assessed. 
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An LPC stress analyst was aware of a reproducible crack occurring at the 
interface between an Inhibited surface an an uninhibited surface. This occurred 
In an Aerojet Mlnuteman motor. The motor had been stress-analyzed in Its Initial 
configuration, but not when a fillet was formed from burning at the Interface: 

fillet 


'-uninhibited V 


uninhibited V- Inhibited 


c 






uninhibited 




inhibited 


The problem manifested itself as a motor failure some time into the bum, and there- 
fore, warranted detailed scudy. (Such Is not the case with the 2.7S inch motor In 
question, which to all Intents and purposes, performs as required.) Under separate 
contract, UTC performed a stress analysis of the Mlnuteman Motor taking the fillet 
Into account and the problem was subsequently Identified and resolved. 


A series of curves was then constructed for an Idealized cylindrical 
2.7S Inch motor and are shown in Figure 6-4. Die first curve shown for the absence 
of port Inhibition, or complete Inhibitor failure on the poet. The second is for a 
partial Inhibition, or a part failure. The third Is for port inhibition per the Aerojet 
design without failure. Die fourth curve takes ihe Aerojet design without Inhibitor 
failure, and adds a diagonally circumferential crack beginning at the fillet end 
propagating at a 45° angle to the wall. In each case, the Is adjusted for the 
measured throat erosion by assuming a constant rate with web burned. 


Each Kn curve Is normalized to the Initial Kjj of the Aerojet design 
without a crack or Inhibitor failure. It Is verified that the port inhibition does 
not significantly alter the progressivlty: the design Inhibition is about 10% more 
progressive than no port inhibition. However, the assumed Ideal crack does 
produce a near-neutrallty. Moreover, the Initial Kjg was determined to be consis- 
tent with the Initial pressure of the motor assuming the predicted erosive burning 
to be valid and not a super erosive burning. 

The ballistics prediction for the Aerojet motor is qualitatively compared 
to the pressure trace furnished by AFRPL In Figure 6-5 , The AFRPL motor has the 
3 ame grain design as the Aerojet motor, but differs from the Aerojet motor as to 
propellant, nozzle configuration, and complete lack of grain inhibition. The follow- 
ing significant aspects may be noted: The AFRPL motor show3 progressivlty: It 

•Ratio of propellant bum area to nozzle throat area. 



Figure 6-4 Diagnostic K N Curves For 2. 75 Motor 
Including Measured Throat Erosion 



Fractional Web Burned 



Figure 6-5 Comparloan Of Predicted Pressure Trace For Aerojet 
2.75 Motor With Actual Trace For AFRPL 2.75 Motor 
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was verified that its progressivity is consistent with its change In Km ; the sharpness 
of its tail-off is closer to the predicted Aerojet motor than the actual Aerojet motor. 
The AFRPL motor shows less progressivity than the predicted Aerojet motor because 
of its complete lack of grain inhibition. It is apparent that the AFRPL motor does 
not exhibit any super erosive burning with high G-dependence . and does not 
exhibit any anomaly. 

For all of the foregoing reasons, it is concluded that the Aerojet 2.75 
inch motor exhibits some reproducible anomaly, that this anomaly is probably a 
crack emanating from the inhibitor Interface, and therefore, it is impossible to 
predict the observed ballistics. On the other hand, the ballistics of the AFRPL 2.75 
inch motor would probably have been successfully predicted were that the test case. 

Since a satisfactory ballistics solution for the 2.75“ motor could not be 
obtained, quantitative comparisons with the measured pressure and thrust versus 
time, and total Impulse, ate not indicative of the accuracy of the program, and are 
not given. A valid comparison with the delivered specific impulse was obtained, 
but in a somewhat non-aprlori manner. Using one of the many operational modes 
of the SPP program, the measured average chamber pressure and erosion rates were 
input, together with an estimate for the average motor L*. The Theoretical Per- 
formance (ODE), Two-Dimensional Two-Phase (TD2P), and Boundary Layer (TBL), 
Modules were then executed. At the high pressure for this motor the execution 
of the Kinetics Module (ODK) is not warranted as the loss is quite small, but the 
amount of computer time that would be required is relatively large. The simplified 
kinetics loss correlation should be adequate at pressures of 1000 psl, or more 
for conventional aluminized propellants. 

Preceding In the aforementioned fashion, the predicted delivered I was 
found to be Just 0.3% above the measured value for this motor. The simplified 
performance prediction correlations yielded a value 0.93% less than the measured 
value for this case. Both of these results are quite satisfactory. 

6.4 Sealed Envelope Prediction; C-4 Third Stage 

This prediction was a special challenge because it was a blind prediction 
of a motor containing a complex grain geometry and a high-exponent propellant. 

As would be expected for such a case, the complete capability of the program was 
utilized In the calculation. Th" specific Impulse prediction for this case was also 
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excellent. The theoretically determined delivered specific Impulse was only 0.4% 
higher than the reported motor firing value, while the simplified performance pre- 
diction correlations yielded a value of 0.6% on the low side. , ; 

The predicted pressure trace Is compared to the actual motor data In Figure • 

6-6. The result Is similar to the Extended Delta result In that pressure is under- . j 

predicted during the Initial portion of the firing, and overpredicted at later times. 

Because the pressure exponent Is high, the extent of the pressure departures Is 
greater than In the Delta. The grain geometry calculations appear to be better 
here In that the pressure peak and dip are predicted to occur at a later time, which 
Is consistent with the predicted slower Initial burning. The compensating errors 
produce a near-correct prediction of action time, and the tall-off Is well-predic- 
ted. 

Action time average pressure was underpredtcted by 0.7%, the best result 
of all test cases. Burning rate was underpredicted by 1.3%. The throat erosion 
was underpredicted by 41% of the average throat erosion; fortunately, the mag- 
nitude of the throat erosion was small and hence, was negligible In performance 
prediction. It Is thought that this error In erosion rate stems from extrapolation 
of the seml-emplrlcal throat erosion constants to a new high temperature regime, 
and should be subject to adjustment In future work. Action time average thrust 
was underpredtcted by 1.3%, Total Impulse was underpredicted by 2. !%• this 
resulted from an error in computed propellant weight, on the low side. The ln.’uts 
for the grain geometry had to be prepared from a small, relatively crude grain 
drawing, and strand burning rate was provided to us at only one condition. In 
view of th.s, the accuracy of these ballistics predictions Is quite good. 

A summary of the results for this case may be found in Table 6-1. A break- 
down of the individual losses, and a comparison between the theoretical and sim- 
plified predictions Is presented In Table 6-2. 

6.S Summary of Results 

The results of the four previously described comparisons with motor firing 
data are summarized In Table 6-1. E t for the anomalous pressure behavior 
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TABLE 6-1 SUMMARY OF RESULTS 





Motor used In developing or evaluating the correlation 



Figure 6-6 Predicted And Actual Pressure Traces For Blind Prediction Test Case 
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observed in the 2.75" motor (discussed in Section 6.3) all of the predicted results 
were well within the accuracy goals for the program. 
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A breakdown of the individual losses for the two cases m which the com- 
plete methodology was executed (Extended Delta and C-4 Third Stage) is pre- 
sented in Table 6-2. The results ot both the theoretical and simplified methods 
are indicated. In this way the relative magnitudes of the different losses car. be 
assessed. The erosion loss, normally tncluded in ( see equation 3-2) has 

been listed separately for the same reason. In order to facilitate comparison be- 
tween the theoretical and simplified 2D-2phase losses, the simplified 2D and 2 
phase losses have been combined (the theoretical loss cannot be separated into 
Us components). 

As a result of the nozzle throat geometry restrictions of the TD2P module, 
the actual nozzle geometry could not be analyzed for the three high performance 
motors: Extended Delta, SRAM, and C4 Third Stage. The transonic flow solution 
is the TD2P module is limited to a simple circular arc th.oat geometry and the 
ratio of the throat radius of curvature to the throat radius must be greater than 
1.5. The two-phase flow loss is usually the largest single loss for highly al- 
uminized propellants, and the magnitude of this loss is sensitive to the throat 
geometry. Therefore, the inability of the TD2P module to treat the actual rozzle 
geometry is cause for some concern. An effort to remove this limitation would be 
warranted. In the mean time, the magnitude of the error resulting from the use 
of a rationally selected approximate nozzle geometry should be ascertained. 
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BroaV-down Of Individual Losses 







































7. CONCLUSIONS AND RECOMMENDATIONS 
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7 .1 Conclusions 

The goal o £ this effort was to develop a computerized performance predic- 
tion methodology that was capable of predicting delivered specific Impulse to +2%, 
and delivered thrust and total Impulse to +5%, for a wide range of solid rocket 
motor geometries , propellmt, and operating conditions. 

The Solid Rocket Performance Program (SPP) described herein, and in the 
other two volumes of this report, was designed to achieve the goal. Six existing 
•computer programs we re integrated with a central logic to form a strong analytical 
foundation for the program. To supplement the theory, where necessary, and to 
increase the flexibility of the program, a number of existing and newly developed 
semi-empirical correlations were Incorporated into the program. 

The program is modular in structure, and allows the user to exercise the 
various modules individually, all together, or in any predetermined combination, 
as required. For those cases where speed and simplicity are paramount (e.g. large 
parametric studies} the analytical portions of the program may be bypassed completely 
in favor of the performance loss correlations which have been incorporated into 
the program. Alternatively, one, or more of these correlations may be employed 
In lieu of the corresponding analytical calculation when the particular performance 
loss ii elthe' known to be small, or adequately characterized by the correlation. 

Ballistics and performance predictions, obtained with the SPP Program, were 
compare with firing data for four widely varying full scale motors. Except for the 
anomalous . .essure-tlme behavior of one of the motors Can explanation for which is 
proffered), allot the predicted results were wi.iln the accuracy goals originally 
set for the program. In each case the difference between the predicted and measured 
delivered specific impulse was less tnan one percent. This Included one case which 
was predicted to within 0.4% on a sealed envelope basi3. The simplified performance 
loss correlations also gave good results for each case. Conclusive statements re- 
garding thr, accuracy and range of validity of the SPP program cannot be made until 
additional verification efforts are conducted. The following tentative conclusions 
have been drawn: 
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»The accuracy of burning rate prediction Is comparable to, and will be limited 
w by, the error In measuring strand burning rates. 

Action time average thrust and pressure appear to be predictable to within 
•sit; compensating errors In the various parameters entering Into the compu- 
tations of pressure and thrust being helpful In this regard. 

^Delivered specific Impulse appears to be predictable to within 2 %, and pro- 
bably to within lit In most cases which do not Involve unusual circumstances. 

0The simplified method of computing specific Impulse efficiency Is adequate 
for Its Intended purpose. 

As a result of the limitations of the present analysis (see Section 7.2) the 
accuracy of the calculated results may be Influenced by compensating errors. For 
Instance, the geometric limitation of the TD2P analysis that requires the actual 
nozzle shape to be approximated In many cases (Including 3 of the present cases), 
the possible errors In particle size distribution, the neglect of throat roughness, 
etc., all Impact the calculated two-phase flow loss. No attempt to quantify the 
magnitude of the possible errors resulting from these and other assumptions was 
made during the course of this Investigation. Further verification efforts, para- 
metric studies and future experimental and analytical investigations will undoubt- 
edly be Illuminating In regards to the adequacy of the individual assumptions and 
the program as a whole. 

The SPP program treats the complete solid rochet motor performance problem, 
I.e., grain geometry, ballistics and efficiency, within the context of a single, 
automated, entity. As a result, accurate analytical predictions of solid rocket 
motor performance can be obtained In a single computer run. It Is hoped that the 
convenience of the present approach will help stimulate more widespread applica- 
tion of theoretical techniques In solid rocket motor development programs. 


Recommendations 


The SPP program has many virtues; it also has its limitations. Recognition 
of these limitations serves a dual purpose: first, to allow the program, as it exists, 
to be more fruitfully employed; secondly, to yield a better perspective on the nature 
of the Improvements and/or extensions of the program that should be considered. 

The following Is a list of the more significant limitations of the SPP program. 


7-2 





1. The TD2P currently terminates If particle Impingement occurs In the 
nozzle. 

2. The two-dimensional two-phase transonic analysis Is inadequate In a 
number of respects: the analysis will fall for normalized radii of cur- 
vature <4.S and for Inlet angles the wall geometry Is limited to 
a conical Inlet and single radius of curvature throat. 

3. The TD2P program also: neglects real gas effects; cannot handle 
shock waves; does not explicitly allow no particle cases to be computed 
(the program can be “fooled" Into doing so), prohibits gas phase con- 
densation and particle growth In the nozzle. 

4. The present boundary layer analysis (TB1) Is not deemed to be accurate 
enough when boundary layer losses are relatively largo. 

5. A method for determining nozzle wall temperature has not been Integra- 
ted Into the methodology. 

6. The nozzle erosion correlations In the program are limited, and not as 
accurate as would be desired. 

7. Throat roughness, and erosion Induced discontinuities are not addressed. 

8. The analysis Is Inadequate In transient dominated cases. 

9. The grain design analysis is Inefficient for "simpler" grain geometries. 

10. The effects of acceleration on burn rate, and nozzle flow separation 
have not been considered. 

11. Three dimensional flow effects have not been considered. 

The State-of-the-ort In solid rocket performance prediction was recently 

(3) 

reviewed at the last JANNAF Performance Standardization Working Group meeting . 
As 3 result of this meeting seven deficiency areas were Identified. These areas 
Included some, or all, of the Items listed under numbers 2, 3, 7 and 11, above. 

In addition, they Included more adequate particle size determination, and better 
bum rrtc modeling. So-called "deficiency statements" were recently prepared by 
several of the comml.iee members, In which the nature of these limitations, and 
methods for alleviating them are ...plored in detail. Based on these statements, 
and our own predilections, it Is recommended that any future efforts to modify the 
SPP program should begin with items 1, 2 and 7, as listed previously. Item 1 Is a 
relatively simple matter to take care of, Items 2 and 7 would require significant 
effort, but are entirely possible to overcome with currently available numerical 
techniques. 
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Most of the other Items In this list represent limitations which are encount- 
ered less frequently (e.g., large boundary layer losses, nozzle separation, tran- 
sient dominated motors, etc.), or ones which would require Inordinate effort to 
overcome (e.g., 3-D effects). Some of these limitations can currently be over- 
come through the use of existing analyses on an uncoupled (to SPP) basts, (e.g., 
boundary layer and tl nrat erosion). Which of these remaining limitations will be 
addressed will most likely be a function of the predilections of future users, and 
the requirements of future motors. 

In passing, one other recommendation should be made. Ihe existing know- 
ledge regarding the kinetics of typical metallzed solid propellant exhaust species 
is, as pointed out In Appendix A, meager, at best. In most cases this Is of little 
consequence In conventional aluminized propellants, as the kinetics loss tends to 
be small. At low pressure, however, the kinetics loss can be significant, and un- 
less additional work towards Identifying the controlling reactions, and their rates, 
Is forthcoming, substantial predictive errors can be expected In such cases. 
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Appendix A - Reaction Rate Screening 

The reaction rate screening done under this contract should by no means 
be considered complete either from an absolute, or from a relative basis. What 
was attempted was to develop a preliminary reaction rate set which would allow 
for the accurate prediction of the losses In performance of a solid rocket motor 
which are due to finite rate gas phase chemical kinetics. 

The criteria by which the reaction rate screening was done was to select 
the minimum number of species and reactions which would produce an acceptably 
accurate value of I sp efficiency ils efficiency Is defined as 

sp KIN 

"KIN ~ T 


where I = I as calculated In the presence of finite rate chemistry 
sp KIN air calculated by the ODK module. 

I =1 as calculated by the restricted equilibrium option of 
sp RE tffi ODE module 

However, at the time the reaction set screening was done, the restricted 
equilibrium option for the ODE module was not completed. This was unfortunate 
since the difference In I gp as calculated by the two options In the ODE module 
can exceed I % and thus Influence the number of species and reactions which could 
dropped from reaction set. 

The basic steps in developing »ha master reaction set for alumlnumlzed 
propellants and the screening of that reaction set are outlined below. 

1. The reaction set of Reference A1 (AP-PBAN with no Al) was merged 

with those reactions in Reference A2 (which Included aluminized species). 
This process yielded a Initial master reaction rate set of 31 species 
and 126 reactions. 

2. The Initial reaction rate set was compared with the results of chemical 
equilibrium calculations and was found to be deficient. 

3. A short literature search was conducted for measured rates concerning 
aluminized species. No measured rate data was found. 

4. Estimates were made for rates Involving 6 more species and 30 more 
reactions using the procedure reported In Reference Al (the same pro- 
cedure was used In Reference A2 to estimate the reaction rates of 
aluminized species). 
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5. Using the Generalized Kinetics Analysis Program (GKAP) of Reference 
A3, the total master reaction set was screened. This computer program 
SCREENED the reaction in an absolute sense. That is, only reactions 
which had negligible influence on the fluid flow or rate of production 
of a chemical species everywhere in the nozzle expansion were drop- 
ped. The total reaction set was reduced only minimally by this first 
step in the screening process. 

6. The next step in screening procedure was to reduce the number of species 
and hence the number of reactions by making an ODE run omitting species 
thought do dbly to be unimportant and comparing the results to an ODK 
calculation using the reduced reaction s ' . This procedure reduced 

the reaction set to 23 species and 73 reactions. 

While the above screening procedure reduced both the number of species 
and reactions considered in performance calculations by the SPP code, it is by no 
means considered adequate. Too many species and reactions were Included in 
this set, which are due to the comparison of two Inconsistent physical models of 
the expansion of a chemically reacting gas through a rocket motor nozzle. Thus, 
the screening done to date should not be considered accurate until measured data are 
taken on the reactions considered in this set and/or performance losses due to 
finite rate gas phase chemical kinetics can be quantified. 

The final reaction set selected for conventional aluminized propellants is 
shown in Table 2-13 of Volume III. The following notes concerning the comment 
field of this reaction set are Intended to help the reader understand where these 
reactions and rate data came from. 

• The first part of each section (third body and binary exchange reactions) 
is from Reference Al. For example, BAULCHO 968)1.1, refers to the mea- 
surements made by Baulch in 1968 and reported in LEEDS 1 report, 

•The second part Is from Reference A2 and refers to the reaction number 
built Into that computer program. 

•The third part, always suffixed by the word ESTIMATE, is from reactions 
and rates estimated by the authors of this document. 
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